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ABSTRACT 


This  report  presents  a  summary  of  the  research  goals  and  accomplishments 
of  a  program  concerned  with  the  dynamic  magneton] ast.ic  properties  of  materials 
containing  rare-earth  elements.  Results  are  presented  for  two  main  classes  of 
phenomena:  first,  the  effects  of  magnctoclastic  interactions  on  the  propagation 
of  elastice  waves  through  rare-earth  materials,  and  second,  the  utilization  of 
the  magnetoelastic  interaction  for  the  generation  of  high-frequency  ultrasonic 
elastic  waves.  Results  were  obtained  for  both  the  pure  rare-earth  elements 
(Gadolinium,  Terbium,  Dysjjrosium,  Holmium,  and  Erbium)  and  for  intermetallic 
compounds  such  as  the  terbium-  iron  cubic  Lavcs-phase  material,  TbFe^.  Attempts 
were  made,  where  possible,  to  interpret  the  experimental  results  in  terms 
of  theoretical  models  for  the  fundamental  magnetoelastic  interactions.  Most  of 
the  work  described  here  was  in  a  preliminary  state  at  the  time  this  contract  was 
terminated,  and  so  the  full  analysis  of  the  work  and  the  ways  in  which  it 
meets  the  goals  of  the  contract  cannot  be  fully  made  at  this  time. 
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I.  INTRODUCTION 

Tliis  report  describes  the  results  of  a  research  program  on  the  magnetic, 
elastic,  and  magnetoelastic  properties  of  1  are-earth  materials,  including  both 
the  pure  rare-earth  elements  and  rare-earth-iron-group  intormetall ic  compounds. 
This  program  has  been  sponsored  by  the  Advanced  Research  Projects  Agency  of  the 
U.  S.  Department  of  Defense,  and  it  has  been  carried  out  under  Contract  Nc . 
DAAH01-71-C-0259  (1  November  1970  -  3]  October  1971)  and  Contract  No.  DAAH01- 
72-C-0285  (4  November  1971  -  31  December  1973),  monitored  by  the  U.  S.  Army 
Missile  Command.  In  this  introductory  section  of  the  report,  the  original 
motivation  for  this  program  and  its  specific  goals  are  reviewed,  and  a  summary 
of  the  accompli sliments  attained  during  the  course  of  the  program  is  presented. 

A .  Motivation  for  Research  on  the  Magr.c1, oclast.ic  Properties  of  Rare-Earth 
Materials . 


As  the  fundamental  properties  of  rare-earth  materials  first  began  to  be 
revealed  through  many  different  types  of  experimental  research,  it  was  quickly 
recognized  that  their  magnetic  and  r. agnctoelastic  properties  are  particularly 
interesting,  with  a  number  of  potentially  important  technological  applications. 

Most  of  the  known  magnetic  and  magn_toolastic  properties  of  the  pure  rare-earth 

« 

elements  and  Pertain  rare-earth  alloys  are  described  in  detail  in  a  recent 

book  edited  by  Elliott'*’,  and  a  compilation  of  many  fundamental  magnetic  and 

magnetoelastic  properties  of  a  number  of  ’-are-earth  materials  is  aiven  by 
2 

Tebble  and  Craik  .  Both  of  these  works  provide  extensive  bibliographies  of  the 
original  literature  in  the  area  of  rare-earth  magnetism. 

One  of  the  most  unusual  maqnetic  properties  of  the  heavy  rare-earth  elements 
((id,  Tb,  Dy ,  Ho,  Fr,  and  Tm  is  the  existence  of  several  types  of  periodic 
magnetic  ordering,  including  both  ferromagnetic  and  anti  ferromagnetic  phases. 
Examples  of  this  periodic  ordering  are  shown  schematically  in  Fig.  1.  All  of 
these  elements  possess  a  hexagonal  close-packed  crystal  structure  (hep)  which 
can  be  considered  to  be  composed  of  basil-plane  layers  jf  trivalent  ions  perpen¬ 
dicular  to  the  hexagonal -symmetry  axis  (c-axis) .  In  all  the  observed  periodic 
ordering  structures,  all  the  ions  in  a  given  basal-plane  layer  are  ferromagneticallv 
aligned,  but  the  direction  of  the  magnetic , moment  of  a  basal-plane  layer  changes 


DAAIIO1-72-C-0285  FINAL  REPORT 


2. 


Fig.  1  Magnetic  Ordering  Arrangements  of  the  Rare-Earth  Elements 


m  a  periodic  fashion  from  one  layer  to  the  next.  For  example,  in  the  case  of 
Dy  .n  its  antiferromagnetic  phase  (85  K  <  T  <  179  K) ,  the  ordering  has  been 
described  as  helicoidal,  with  the  moment  of  each  basal-plane  layer  lying  in  the 
basal  plane,  but  with  a  finite  angle  of  rotation  of  the  magnetic  moment  from 
one  layer  to  the  next,  as  indicated  in  Fig.  1. 

The  periodic  magnetic  ordering  in  these  rare-earth  elements  arises  from  the 
nature  of  the  exchange  interaction  between  rare-earth  ions,  an  indirect  interaction 
in  which  the  conduction  electrons  play  an  essential  part.  This  indirect  exchange 
interaction  may  be  envisioned  as  one  in  which  a  rare-earth  ion  interacts,  through 
the  normal  electrostatic  exchange  interaction,  with  the  conduction  electrons 
at  the  Fermi  surface,  leading  to  a  spin  polarization  of  the  sea  of  conduction 
electrons  in  the  vicinity  of  the  rare-earth  ion.  This  spin  polarization  extends 
through  the  crystal  as  a  spin-density  wave,  with  a  wavelength  characteristic  of 
the  electronic  wavelength  for  electrons  at  the  Fermi  energv.  The  conduction  elec¬ 
trons  then  undergo  direct  exchange  interactions  with  neighboring  rare-earth  ions, 
and  it  has  been  shown1  that  the  net  interaction  between  two  rare-earth  ions  is 
such  that  it  can  be  described  by  means  of  a  Heisenberg  exchange  operator.  Because, 
however,  the  electron  polarization  produced  by  the  interaction  with  the  first 
rare-earth  ion  oscillates  spatially  in  a  wavelike  fashion,  the  interaction 
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between  two  neighboring  rare-earth  ions  may  be  either  ferromagnetic  or  antiferro- 
magnetic,  depending  upon  their  spatia1  separation.  It  is  the  competition  between 
anti  ferromagnetic  and  ferromagnetic  exchange  among  all  pairs  of  ions  in  the  crystal 
which  leads  to  periodic  magnetic  ordering.  The  actual  type  of  ordering,  however, 
is  a  result  of  both  the  i xchange  interactions  and  the  magnetic  anisotropy  of  a 
particular  material. 

Although  the  phenomenon  of  periodic  ordering  associated  with  the  pure  rare-eart 
elements  and  many  of  their  alloys  and  comoounds  is  in  itself  quite  interestina, 
it  does  not  constitute  the  principal  point  of  interest  from  the  standpoint  of  the 
potential  technological  applications  of  rare-earth  materials  which  are  the  subject 
of  the  research  program  described  here.  The  periodic  ordering  found  in  many  of 
these  materials  does,  however,  play  an  important  part  in  the  determination  of 
many  other  observable  properties  of  interest,  and  it  must,  of  course,  be  carefully 
taken  into  account  when  it  occurs . 

Most  of  the  properties  of  interest  for  the  present  research  program  are 
connected  with  the  very  large  magnetocrystalline  anisotropy  exhibited  by  most 
magneticallv  ordered  rare-earth  materials  and  with  the  dependence  of  this  aniso¬ 
tropy  on  the  state  of  strain  of  the  crystal  lattice.  The  strain  dependence  of  the 
anisotropy  leads,  of  ccurse,  to  the  phenomenon  of  magnetostriction.  Although 
an  appreciable  fraction  of  the  overall  magnetic  anisotropy  of  rare-earth  magnetic 
materials  arises  from  the  anisotropic  nature  of  the  exchange  interaction  itself, 
the  most  important  contribution  arises  from  the  interaction  of  each  rare-earth 
ions  witli  the -crystalline  electric  field  produced  by  the  neighboring  ions. 

Although  the  crystalline  electric  field  is  small  in  metals,  because  of  the  shielding 
effect  of  the  conduction  electron,,  its  effect  on  the  orientation  of  the  angular 
momentum  of  a  rare-earth  ion  is  quite  large.  The  electrons  responsible  for  the 
magnetic  moment  of  a  rare-earth  ion  lie  in  the  r-sholl,  highly  localized  deep 
within  the  ion  and  shielded  from  the  environment  by  filled  outer  s-  and  p-shells. 
Because  the  f-electrons  retain  very  closely  the  same  character  in  a  crystalline 
environment  as  that  found  in  free  ions,  the  electronic  charge  distribution  of 
the  f-shell  of  a  rare-earth  ion  is  highly  anisotropic.  Consequently,  the  free 
energy  associated  with  each  ion  in  a  crystal  depends  strongly  upon  the  orientation 
of  the  angular-momentum  vector  of  the  ion  with  respect  to  tha  crystallographic 
principal  axes  associated  with  the  crystalline  electric  field.  Thus,  the  magnetic 
properties  of  ’-are-earth  materials  are  nearly  always  highly  anisotropic,  with 
effective  uniaxial  anisotropy  fields  as  large  as  several  hundreds  of  kilogauss. 


A 
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The  extremely  large  anisotropy  of  the  rare-earth  materials  has  led  to  the  develop¬ 
ment  of  a  new  class  of  rare-earth-cobalt  permanent  magnets  with  very  large 
coercive  force  and  high  remanent  magnetization  .  For  the  research  program  described 
hero,  however,  it  is  not  so  much  the  magnitude  of  the  anisotropy  energy  that  is  of 
interest  as  it  is  the  dependence  of  this  anisotropy  energy  upon  the  state  of 
lattice  strain. 

Because  the  magnetic  anisotropy  depends  largolv  upon  the  crystal  field  pro¬ 
duced  by  the  ions  in  the  neighborhood  of  a  given  lattice  site,  and  becau-e 
the  magnitude  and  symmetry  of  this  crystal  field  depend  strongly  upon  the  state 
of  strain  of  the  lattice  with  respect  to  its  normal  unstrained  state,  the  overall 
free  energy,  a  combination  of  the  anisotropy  energy  and  the  elastic  energy,  may 
reach  a  minimum  for  a  finite  value  of  the  lattice  strain  tensor,  leading  to  the 
phenomenon  of  magnetostriction.  In  the  case  of  many  rare-earth  materials,  the 
observed  magnetostrictive  strain  can  be  as  large  as  several  orders  of  magnitude 
greater  than  that  observed  in  most  iron-group  magnetic  materials.  This  extremely 
large  magnetostriction,  manifested  by  magnetostrictive  strain  as  large  as  nearlv 
one  per  cent  in  some  mat  trials  ,  is  the  property  of  rare-earth  materials  which 
has  motivated  most  of  the  research  program  described  here.  The  fact  that  the 
magnetostriction  observed  in  many  rare-earth  materials  is  many  times  larger  than 
that  of  most  other  magnetostrictive  materials  suggests  the  possibility  that 
rare-earth  materials  may  provide  a  new  class  of  materials  suitable  for  applications 
as  magnetostrictive  ultrasonic  transducers  of  high  efficiency.  Such  transducers, 
capable  of  the  generation  of  high-intensity,  high-frequency,  sound  waves  might  be 
expected  to  lead  to  significant  improvement  in  the  development  of  sonar  trans¬ 
ducers  and  in  the  area  of  materials  testing  and  processing.  Most  of  the  work 
described  here  has  been  directed  toward  the  obtaining  of  a  better  understanding 
of  magnetoelastic  effects  in  rare-earth  materials,  particularly  in  the  case  of 
dynamic  maanetoelastic  strain,  a. id  toward  the  development  of  high-frequency 
ultrasonic  transducers  utilizing  rare-earth  materials. 

In  addition  to  the  possible  application  of  the  large  magnetostriction  of 

rarr-earth  materials  in  the  development  of  highly  efficient  ultrasonic  transducers, 

work  on  the  general  magnetoclastic  properties  may  prove  useful  in  providing 

information  of  importance  in  the  development  of  rare-earth  permanent  magnets^ 

and  in  the  development  of  rare-earth-iron-group  magnetic  bubble-domain  memorv  and 
4 

logic  elements  .  Furthermore,  rare-earth  materials  may  find  useful  applications, 
as  a  result  of  some  of  the  work  described  here,  as  acoustic  circuit  elements  in 
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hiqh- frequency  signal-processing  systems  employing  both  bulk  and  surface  acoustic 
waves. 

B.  Objectives  of  the  Research  Program 

In  the  initial  proposal  which  led  to  the  funding  of  this  research  program 
by  the  Advanced  Research  Projects  Agency,  the  principal  objective  was  the  develop¬ 
ment  of  highly  efficient  high-frequency  r agnetostrictive  ultrasonic  transducers 
utilizing  rare-earth  materials.  Although  it  was  realized  that  the  pure  rare-earth 
elements,  which  are  strongly  magnetostrictive  only  below  their  magnetic  transi¬ 
tion  temperatures,  nearly  all  of  which  are  in  the  cryogenic  range,  would  probably 
not  be  useful  in  many  applications  where  the  cryogenic  temperatures  required  for 
their  operation  as  transducers  would  be  unacceptable,  it  was  initially  proposed  to 
concentrate  a  considerable  effort  on  transducers  fabricated  from  these  elements 
in  order  to  learns  as  much  as  possible  concerning  dynamic  magnetoelastic  effects 
at  high  frequencies  in  rare-earth  materials.  The  ultimate  goal,  of  cour  e,  was 
to  find  new  materials  with  large  magnetostriction  at  room  temperature  and  above. 
Such  materials  were  unavailable  at  the  beginning  of  this  program,  but  in  1971 
Clark  and  Belson^  reported  the  observation  of  very  large  magnetostriction  in 
polycrystalline  TbFe2  and  other  rare-earth-iron  interne tallic  compounds.  From 
that  point  in  time,  the  research  carried  out  in  this  program  on  the  development 
of  magnetostrictive  ultrasonic  transducers  was  concentrated  on  these  new 
materials . 

In  additior  to  work  on  the  development  o:  ultrasonic  transducers,  it  was 
apparent  from  the  inception  of  this  research  program  that  it  would  be  necessary 
to  obtain  information  concerning  a  number  of  fundamental  properties  of  rare-earth 
materials  related  to  the  dynamic  magnetoelastic  effects  of  primary  interest.  In 
particular,  the  effects  of  magn  :toc-lastic  interactions  on  the  characteristics  of 
elastic-wave  propagation  in  rare-earth  materials  are  important  in  that  they  are 
directly  related  to  the  problem  of  the  magnetostrictive  generation  of  ultrasonic 
elastic  waves.  Consequently,  a  program  was  initiated  to  obtain  information  on  the 
elastic  constants  of  the  rare-earth  elements,  with  emphasis  on  the  magnetoelastic 
contributions  to  these  elastic  constants,  and  on  elastic-wave  attenuation  at  hiah 
frequencies  in  the  rare-earth  elements.  At  the  same  time,  a  theoretical  analysis 
of  the'problem  of  magnetoelastic  effects  on  elastic-wave  propagation  was  developed 
in  order  to  permit  a  satisfactory  analysis  of  the  experimental  results  in  terms 
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of  established  theoretical  models  of  the  static  magnetic  and  magnetoelastic 
properties  of  rare-earth  materials.  This  theoretical  treatment  of  dynamic  macnetoela 
tic  effects  in  rare-earth  materials  was  also  formu]ated  in  such  a  way  as  to  cerm.it 
the  prediction  of  the  characteristics  of  magnetostrictive  ultrasonic  generation 
in  these  materials,  but  it  was  limited  to  the  case  of  single-crystal  specimens. 

Othei  experimental  programs  wore  developed  as  the  overall  program  progressed, 
for  the  purpose  of  correlating  static  magnetic  and  magnetoelastic  properties  of 
rare-earth  materials  with  their  dynamic  magnetoelastic  properties.  One  of  these 
programs  concerned  the  measurement  of  the  magnetization  of  thin-film  specimens  of 
the  materials  used  for  the  fabrication  of  thin-film  ultrasonic  transducers.  It 
was  discovered  very  early  in  th'  use  of  TbFe2  for  ultrasonic  transducers  t  .at  the 
properties  of  this  material  in  the  form  of  thin  films  are  quite  different  from  the 

^  it 

reported  bulk  properties.  Consequently,  a  sensitive  vibrating-sample  magnetometer 
was  constructed  for  the  purpose  of  measuring  the  magnetization  of  films  deposited 
on  various  substrates  under  a  wide  range  of  deposition  conditions.  Another 
program  concerned  the  measurement  of  the  magnetostriction  of  thin-film  specimens 
and  very  small  single-crystal  specimens  of  rare-earth-iron-group  intermetaliic 
compounds  such  as  TbFe2>  A  system  utilizing  the  method  of  holographic  interfero¬ 
metry  was  developed  for  this  purpose,  because  it  was  believed  that  only  such 
a  method  would  permit  the  measurement  of  the  magnetostriction  of  fragile  and  brittle 
thin  films  of  rare-earth  materials. 

C.  Principal  Accomplishments  of  the  Resear^n  Program 

The  major  accomplishments  of  this  research  program  on  the  magnetoelastic 
properties  of  rare-earth  materials  are  summarized  briefly  here.  A  more  detailed 
description  of  each  area  is  given  in  the  main  body  of  this  report  which  follows. 

1.  Magnetostrictive  Ultrasonic  Generation:  Magnetostrictive  ultrasonic 
generation  was  demonstrated  in  thin-film  transducers  of  several  pure  rare-earth 
elements,  namely  Gd,  Tb,  Ho,  Dy,  Er,  and  in  thin-film  transducers  of  TbFe2  and 
other  rare-earth-iron  intermetaliic  compounds.  In  the  case  of  the  pure  elements, 
particularly  Dy  at  liquid-helium  temperatures,  very  strong  ultrasonic  generation 
3t  frequencies  up  to  1,550  MHz  was  observed,  with  efficiencies  of  generation 
some  30  dB  better  than  that  observed  in  the  case  of  piezoelectric  generation  in 
single-crystal  quartz  under  the  same  experimental  conditions.  As  expected,  however, 
magnetostrictive  ulti asonic  generation  could  not  bo  observed  at  temperatures  above 
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the  magnetic  ordering  temperatures  of  these  elements.  One  exception  occurred, 
however,  in  the  case  of  Tb,  ir  which  ultrasonic  generation  was  observed  at 
temperatures  in  the  paramagnetic  range  (T  >  230  K) .  In  this  case,  however ,  the 
ultrasonic  intensity  was  fairly  small.  Some  ultrasonic  generation  was  observed 
at  frequencies  as  high  as  10  GHz,  but  only  at  very  low  temperatures,  due  to  the 
large  attenuation  of  elastic  waves  of  this  high  frequency  at  temperatures  above 
20-30  K  in  most  materials.  Because  of  the  interest  in  developing  transducers 
of  high  efficiency  capable  of  operation  at  room  temperature  and  above,  work 
on  transducers  fabricated  from  the  pure  rare-earth  elements  was  not  continued 
after  it  became  evident  that  the  rare-carl h-iron  intcrmctallic  compounds  were 
considerably  more  promising.  Nevertheless,  as  described  more  fully  in  the  fol¬ 
lowing  section,  a  gret  amount  of  useful  information  concerning  ultrasonic  generation 
in  the  rare-earth  elements  was  obtained.  Although  it  was  originally  planned  to 
investigate  magnetostrictive  ultrasonic  generation  in  single-crystal  specimens  of 
the  pure  rare-earth  elements,  it  was  not  possible  to  obtain  results  .in  this  area 
during  the  time  span  of  this  program. 

The  most  interesting  results  were  obtained  in  the  case  of  magnetostrictive 

ultrasonic  generation  in  thin- film  transducers  fabricated  from  the  intermetallic 

5 

compound  TbFe^.  As  mentioned  above,  Clark  and  Belson  found  that  this  material 
in  bulk  form  possesses  a  very  large  static  magnetostriction  at  room  temperature 
in  polycrystalline  form.  This  result  is  in  strong  contrast  to  the  case  of  the 
pure  rare-earth  elements,  for  which,  even  in  highly  magnetically  ordered  states 
at  low  temperatures,  large  magnetostriction  is  found  only  in  single-crystal, 
specimens.  Thin-film  transducers  evaporated  onto  single-crystal  quartz  substrates 
were  fabricated  in  this  proaram,  and  it  was  quickly  found  that  strong  magnetostric¬ 
tive  ultrasonic  generation  occurs  at  frequencies  as  high  as  1,550  MHz.  Only 
shear-wave  generation  was  observed,  a  result  accounted  for  by  the  theoretical 
model  presented  in  a  subsequent  section,  but  the  intensity  of  the  ultrasonic  waves 
generated  magnetostrictivcly  was  some  20  dB  greater  than  that  of  piezoelectrically 
generated  longitudinal  waves  in  single— crystal  quartz  and,  perhaps,  30  dB  greater 
than  that  of  piezoelectri call”  generated  shear  waves  at  the  same  temperature  and 
frequency . 

An  unexpected  but  highly  interesting  result  of  the  investigation  of  magneto¬ 
strictive  ultrasonic  generation  in  thin  films  of  TbFe^  was  the  occurrence  of  a 
very  large  hysteresis  in  the  field  dependence  of  the  ultrasonic  intensity.  After 
initial  magnetization  to  saturation  in  an  applied  field  of  approximately  20  kG, 
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the  film  still  exhibited  strong  ultrasonic  generation  when  the  field  was  removed. 

An  examination  of  the  film  magnetization  using  a  vibrating-sample  magnetometer 
revealed  that  the  film  _cts  as  a  permanent  magnet,  with  a  coercive  force  of 
approximately  4  kOe  at  room  temperature  and  a  remanent  magnetization  of 
approximately  80  pier  cent  of  the  saturation  value.  This  thin-film  behavior  is  in 
strong  contrast  to  that  of  bulk  TbFe2,  which  exhibits  fairly  small  hysteresis, 
with  a  correspondingly  small  coercive  force  and  remanence.  Although  this  interesting 
behavior  is  not  understood  at  present,  it  is  perhaps  due  to  the  occurrence  of  some 
degree  of  crystallographic  ordering  in  the  film  during  the  deposition  procedure. 

One  interesting  possible  application  of  this  large  coercive  force  and  remanence 
is  the  development  of  thin-film  permanent  magnets  which  may  be  useful  in  magnetic 
bubble-domain  memory  and  logic  systems  requiring  permanent  bias  fields  over  a  very 
small  volume. 

At  the  time  of  writing  this  report,  work  is  still  in  progress  on  the 
properties  of  thin-film  'ltrasonic  transducers  fabricated  from  TbFe2  and  other 
rare-earth-iron  compounds.  The  effects  of  variation  of  the  deposition  parameters, 
film  thickness,  composition,  and  the  effects  of  post-deposition  annealing  are 
being  systematically  studied.  Because  of  the  totally  unexpected  hysteresis  encounter 
in  the  work  described  above,  an  investigation,  using  the  method  of  holographic 
interferometry  described  in  a  subsequent  section  of  this  report,  of  the  static 
magnetostriction  of  thin  films,  both  free  and  bonded  to  a  substrate,  is  under  way. 
Since  the  magnetization  exhibits  such  unusual  properties  when  TbFe^  is  produced  in 
the  form  of  a  thin  film,  it  is  likely  that  other  magnetic  and  magnetoelastic  pro¬ 
perties  of  this  material  will  also  be  quite  different  from  their  bulk  behavior. 

2 .  Elastic-Wave  Propagation  in  Rare-F.arth  Materials:  As  mentioned  above , 
an  extensive  program  of  the  investigation  of  the  behavior  of  both  the  elastic-wave 
velocities  and  elastic-wave  attenuation  in  rare-earth  materials  was  carried  out 
under  the  overall  research  program.  Studies  were  carried  out  on  single-crystal 
specimens  of  the  pure  rare-earth  elements  Tb,  Dy,  Ho,  and  Er,  and,  in  addition  to 
providing  in  the  cases  of  Tb,  Ho,  and  Er  new  information  on  the  elastic  constants 
of  these  elements  over  wide  temperature  ranges,  this  work  yielded  a  number  of 
interesting  and  useful  results  concerning  magnetoelastic  effects  on  the  behavior 
of  elastic  waves  at  higli  frequencies  in  rare-earth  materials.  In  particular, 
both  linear  and  circular  magnetoacoustic  birefringence  was  observed  in  these 
materials,  in  accordance  with  the  predictions  of  the  theory  described  in  the  fol- 
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lowing  section.  Those  birefringence  results  permit  the  determination  of  certain 
of  the  shear  magnetoelastic  coupling  constants  of  rare-earth  materials  in  a 
manner  not  possible  without  considerable  difficulty  by  other  meins,  and,  more 
important,  they  permit  this  determination  at  frequencies  and  temperatures 
of  interest  for  magnetostrictive  ultrasonic  generation.  It  is  the  shear  magneto- 
elastic  couplings  which  are  of  greatest  importance  for  the  magnetostrictive 
generation  of  ultrasonic  waves  at  high  frequencies.  The  observed  birefringence 
effects  also  load  to  the  possibility  of  the  development  of  acoustic  circuit 
elements  with  non-reciprocal  properties,  such  as  isolators  and  circulators. 

It  was  also  observed  in  some  cases  that  elastic-wave  velocities  could  bo 
varied  by  as  much  as  7-0  per  cent  by  the  application  of  a  magnetic  field  under 
appropriat"1  conditions.  Such  an  effect  may  be  useful  in  the  development  of 
simple  acoustic  variable  delay  linos.  In  other  cases,  strong  field-dependent 
ultrasonic  attenuation  was  observed,  leading  to  the  possibility  that  simple 
acoustic  variable  attenuators  could  be  developed  using  rare-earth  materials. 

Many  of  the  results  obtained  in  this  part  of  the  program  could  be  analyzed 
by  means  of  straight forward  theoretical  models  of  dynamic  magnetoelastic  interac¬ 
tions,  including  one  presented  in  the  following  section  and  a  theory  developed 
by  Southern  and  Goodings  6.  The  results  of  this  theoretical  analysis  indicate 
quite  good  agreement  between  theory  and  experiment  in  the  area  of  dynamic  magneto- 
elastic  interactions  in  rare-earth  materials,  and,  particularly,  they  show  that 
measurements  on  the  propagation  of  elastic  waves  in  such  materials  may  provide 
the  best  way  to  determine  completely  ti  e  quantities  which  characterize  the 
magnetoelastic  interactions  which  can  occur. 

3.  Holographic  Measurement  of  Magnetostriction :  A  system  for  the 
measurement  of  stai.ic  magnetostriction  in  rare-earth  materials  employing  the 
method  of  real-time  holographic  interferometry  was  developed  to  permit 
measurements  in  specimens  of  rare-earth  materials  whose  configuration  would  not 
permit  the  use  of  more  conventional  methods  for  the  measurement  of  magnetostriction. 
Such  configurations  include  thin-film  specimens  whose  magnetostriction  may  deviate, 
as  discussed  above,  strongly  from  that  of  the  same  material  in  bulk  form,  and  they 
include  the  very  small  single-crystal  specimens  which  can  be  produced  from  rare- 
earth-iron  compounds .  After  the  confirmation  that  the  holographic  method  does 
-’ndeed  give  results  in  good  agreement  with  those  obtained  *  v  conventional  methods 
(for  specimens  which  permit  the  use  of,  for  exam1  le,  strain-gauge  methods),  an 
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attempt  was  made  to  apply  this  method  to  the  case  of  thin  films  of  TbFe2  and 
to  the  case  of  a  small  single  crystal  of  ErFe^  which  was  grown  for  this  ourpose. 
Because  of  the  fragility  of  the  thin  films,  it  has  not  yet  been  possible  to 
obtain  results  in  this  area.  However,  an  extensive  set  of  results  on  the  static 
magnetostriction  of  the  si nc le-crystal  specimen  were  successfully  obtained.  When 
these  results  showed  considerable  disagreement  with  those  reported  (and  confirmed 
in  this  work)  for  polycrystalline  ErFe2,  the  crystal  was  examined  by  the  method 
of  X-ray  fluorescence,  using  an  electron  microprobe  for  excitation,  in  order  to 
determine  its  composition.  Despite  the  fact  that  the  crystal  had  been  prepared 
from  material  for  which  the  stated  purity  of  the  erbium  component  was  99  per  cent, 
according  to  the  commercial  supplier,  it  was  found  that  the  crystal  was,  in  fact, 
a  mixture  of  ErFe^  and  TbFe2>  with  rouqhly  comparable  amounts  of  Er  and  Tb.  Thus, 
the  crystal  should  not  have  been  expected  to  give  results  agreeing  with  what  was 
approximately  expected  for  pure  ErFe2>  Numerous  cross  checks  did,  however,  confirm 
that  the  magnetostriction  was  measured  accurately,  and  the  result  of  the  investi¬ 
gation  indicates  that  this  TJ"JxFr^_yFe2  crystal  exhibits  room  temperature  magneto¬ 
striction  somewhat  large  than  had  been  reported  by  Clark  and  Belson^  for  poly¬ 
crystalline  TbFe^.  The  analysis  of  the  results  is  in  progress,  and  an  attempt 
to  grow  a  second  } rFe2  crystal  using  material  of  known  high  purity  is  under  way. 
Work  is  also  continuing  on  the  effort  to  measure  the  static  magnetostriction  of 
thin-film  Tbre^  in  order  to  perhaps  obtain  a  better  understanding  of  the 
permanent  magnet-like  behavior  of  these  films. 

Theoretical  Analysis  of  Dynamic  Maqnetoelastic  Effects:  When  this 

program  began,  very  little  attention  had  bc°n  given  to  the  theory  of  dyn.imic 

magnetoelastic  effects  in  rare-earth  materials.  A  general  theory  of  static 

magnetoelastic  effects  developed  by  Callen  and  Callen7  was  at  that  time  well 

established,  and  its  predictions  agreed  quite  well  with  many  experimentally 

determined  magnetoelastic  properties  of  rare-earth  materials,  but  the  only 

attempts  to  explain  dynamic  magnetoelastic  effects  in  these  materials  up  to  that 
•  . .  8,9  .  , 

point  in  time.  had  largely  been  nased  upon  earlier  theoretical  models  developed 
for  application  to  iron-group  materials.  It  was,  thus,  recognized  at  the  outset 
of  this  program  that  a  new  approach  was  needed  in  this  area,  extendinq  the  theory 
of  Callen  and  Callen  to  include  dynamic  magnetoelastic  effects  and  adding  to  the 
established  theory  where  necessary.  Other  researchers  began  to  attack  the  problem 
at  the  same  time.  Freyne  developed  a  thermodynamic  approach  to  explain  certain 
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aspects  of  the  anomalous  temperature  de,  andence  of  the  elastic  constants  of  Gd, 
and  southern  and  Goodings  '12  developed  a  somewhat  more  general  theory  of  the 
magnetoelastic  contributions  to  the  elastic  constant  of  rare-earth  materials. 

The  work  of  Southern  and  Goodings  is  particularly  important,  since  these  authors 
took  account  of  the  intrinsic  nonlinear  nature  of  elastic  strain  and  also  deve- 
loped  a  rotational ly  invariant  model  of  the  magnetoolasf ic  interaction,  including 
terms  which  are  important  in  the  case  of  dynamic  elastic  strain  but  which  had  been 
neglected  in  earlier  treatments  of  static  mannetoelastic  effects. 

In  the  present  work,  as  described  in  the  folio-wing  section,  a  theory  of  dynamic 
magneto  elastic  interactions  in  rare-earth  materials  has  been  developed  in  a  way 
which  combines  parts  of  the  work  of  Freyno10  and  that  of  Southern  and  Goodings11'12 
This  theoretical  treatment  has  been  applied  to  explain  satisfactorily  certain  of 
the  experimental  results  reported  here,  and  it  is  shown  in  the  following  section 
how  it  can  be  extended  to  apply  to  many  more  magnetoelastic  phenomena. 

D.  Summary  of  thn_Progrnn  and  Its  Overall  . 

In  some  ways,  the  research  program  discussed  lure  accomplished  its  overall 
objectives  satisfactorily,  but  there  were  also  ways  in  which  it  fell  short  of 
the  original  goals.  A  number  of  interesting  new  experimental  results  on  the 
dynamic  magnetoelastic  properties  of  rare-earth  materials  were  obtained  and 
reported  in  the  literature  and  at  technical  conferences,  and,  in  particular,  it 
was  shown  that  the  intermetallic  compound  TbFe,  is,  as  expected,  a  very  promising 
material  for  the  development  of  room-temperature  magnctostrictive  ultrasonic 
transducers.  The  time  available  for  this  program  did  not,  however,  oermit  the 
complet-on  of  all  the  individual  research  projects  of  interest.  Much  work  remains 
to  be  cone  on  the  complete  characterization  of  TbFc.,  and  other  rare-earth-iron 
intermetallic  compounds  as  magnctostrictive  ultrasonic  transducers.  For  example 
time  did  not  permit  the  investigation  of  the  effects  of  varying  the  thickness  of 
the  transducer  films,  of  deposition  in  an  applied  magnetic  field,  or  of  post-depo¬ 
sition  annealing  of  the  films,  nor  was  it  possible  to  Investigate  the  properties  of 
a  wide  range  of  film  compositions. 

Further  work  on  the  theory  of  dynamic  magneton lastic  interactions,  including 
an  attempt  to  extend  it  to  cover  very  high  frequencies,  is  also  necessary.  There 
also  remains  a  large  body  of  experimental  data  for  which  there  has  not  been  time 
to  apply  the  theory  to  explain  the  quantitative  results.  It  is  believed  that 
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further  extensive  ciilculiilioiis  dloiuj  these  lines  would  permit.  the  strengthening 
of  the  theory  ami  an  overall  bet  ter  under r;l  andi ng  of  dynamic  maenotoel antic  effects 
in  rare-earth  materials. 

The  A.R.I’.A.  .support  of  this  program  also  ended  at  a  time  when  the  holo¬ 
graphic  system  for  tin:  measurement  of  magnetostrict ion  in  thin  films  and  very 
small  single  crystals  was  just  beginning  to  yield  very  interesting  results. 

It  is  felt  that  this  method  will  very  usefully  supplement  ofchei  methods  for  the 
measurement  of  static  magnet  oelast.ic  effects  when  it  reaches  its  full  potential. 

In  all,  it  can  be  said  that  this  program  produced  a  number  of  interesting 
and  useful  now  results,  but  that  it  ended  at  a  time  when  the  most  interest ing 
stage  of  the  program  had  been  reached.  It  is  hoped,  of  course,  that  many  of  the 
individual  projects  described  above  and,  in  more  detail,  in  what  follows,  can 
he  continued  with  a  new  source  of  funding. 


II.  Til  KORY  OF  DYNAMIC:  MAGNKTOKUUVriC 


INTKlvACTlONf'  IN  RA10 ’-KAR’i’l!  MAVKRTAbS 


Most  recent  attempts  to  develop  a  satisfactory  theoretical  model  of 
dynamic  magnet  ''elastic  effects  in  tar  ’-earth  materials,  have  been  based  upon 
the  ext  ensroi  of  t  no  original  theory  ot  static  magii<  toelastic  interactions  in 
tlu.  i  ai  e-ear  Lh  elements  ol  fallen  anil  fallen  to  include  tin  etlects  of  time— de¬ 
pendent  lattice  s.t  rain.  Freyne  emj'loyed  a  quasi  static  thermodynamic  approach 
to  calculate  t  he  e t  t  ec t  s  of  tin*  m.ig ne t  oe  1  a?d  i c  i  n t  e rae t  i cri ;s  oi.  *  he  elastic 
constants  of  gadolinium,  achieving  quite  sat ist artory  results.  His  approach  is 
quite  similar  to  that  employed  here,  although  the  present  treatment,  described  m 


detail  below  can  be  extended  to  cover  a  wider  range  of  phenomena.  Southern  and 


11  ,12 

Goodings  '  ‘  employed  a  quantum- f i eld-t boor et i c  treatment  to  consider  the  same 


general  class  of  problems,  obtaining  in  some  rases  result}!  essentially  in  agreement 
with  those  of  Freyne.  The  most  interesting  new  contribution  of  Southern  and  Good¬ 
ings,  however,  was  the  introduction  into  the  Hamiltonian  which  serves  as  a  model 
for  the  magnet oel ast i c  int enactions  two  new  features.  First,  they  took  explicit 
account  of  the  tact  that  1 inite  s_.  tins  are  quadratic  in  the  components  of  the 
tensor,  gradient  ol  the  lattice  di:  'placement ,  and,  second,  they  t  ook  explicit 
account  of  terms  in  the  Hamiltonian  needed  in  order  that  it  bo  rotational ly 
invariant.  In  earlier  treatments  ot  static  magnet oelast.ic  t  oets,  in  which  a 
homogeneous ,  uniform,  strain  tensor  was  assumed,  these  additional  terms  do  not 
appear.  They  become  important,  however,  when  lattice  waves  of  t inite  wavelength 
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arc  considered,  and,  according  to  Southern  and  Goodings11,  these  additional 
terms  lead  to  interesting  observable  effects. 

The  following  treatment  makes  use  of  a  number  of  ideas  originated  by  both 
Freyne10  and  by  So  hern  and  Goodings11,  but  it  is  formulated  in  such  a  way  that 
it  can  in  many  cases  be  applied  to  the  calculation  of  specific  effects  with 
less  difficulty,  and  it  is  also  useful  in  the  calculation  of  magnetostrictive 
ultrasonic  generation,  a  topic  not  considered  by  either  Freyne  or  Southern  and 
Goodings.  The  treatment  presented  here  makes  use  of  a  number  of  simplifying 
approximations  for  the  purpose  of  facilitating  specific  calculations,  but  nany 
of  those  simplifications  can  undoubtedly  be  eliminated  in  cases  where  the  added 
difficulty  of  making  calculations  appears  warranted.  The  model  presented  hero  is 
not,  therefore,  intended  to  bo  a  definitive  treatment  of  the  problem  of  dynamic 
magnetoelastic  effects  in  rare-ear.th  materials;  rather,  it  is  a  model  which  permits 
semiquantitativo  calculations  which  are  in  some  cases  in  quite  good  agreement 
with  experimental  results,  and  it  is  a  model  which  can  bo  improved  when  necessary. 

A.  Description  of  the  Model  for  Dynamic  Hagnetoelastic  Interactions 

There  is  very  lit'  ^  difference  between  the  basic  model  of  a  system  of 

rare-earth  ions  in  a  crystal  lattice  necessary  for  the  description  of  dynamic 

effects  and  that  which  has  been  used  by  many  authors  for  the  description  of 

static  magnetoelastic  effects.  Much  of  the  earlier  work  on  static  magnetic  and 

magnetoelastic  properties  of  rare-earth  materials  is  summarized  in  the  recent 

book  edited  by  Elliott1,  and  a  more  detailed  treatment  of  maonetoclastic  interac- 

7 

tions  is  given  by  Callen  and  Callon  .  The  principal  difference  encountered  in 
setting  up  a  model  Hamiltonian  for  the  purpose  of  describing  dynamic  magretoelastic 
effects  is  the  necessity  of  allowing  for  strain  waves  of  finite  wavelength  and 
for  spin  waves  of  finite  wavelength.  Most  of  the  effects,  however,  to  which  the 
present  treatment  has  been  applied  involve  strain  waves  of  wavelengch  much  greater 
than  the  dimensions  of  a  unit  crystallographic  cell,  and,  for  the  purpose  of 
simplifying  calculations,  it  has  been  assumed  in  most  of  what  follows  that  any  dyna¬ 
mic  motion  of  the  magnetic  moments  of  the  rare-earth  ions  can  be  described  by 
means  of  a  uniform-precession  mode  of  motion.  In  effect,  this  approximation 
ignores  the  spin-wave  dispersion  at  long  wavelengths,  and  it  is  probably  very  well 
justified,  since  most  of  the  effects  to  which  the  present  treatment  has  been  applied 
involve  elastic  waves  of  wavelength  of  the  order  of  an  appreciable  fraction  of 
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a  millimetei  Thus,  the  principal  additions  to  the  usual  static  model  Hamiltonian 
are  simply  the  introduction  of  the  lattice  kinetic  energy  associated  with  strain 
waves  of  finite  frequency  and  che  introduction  of  the  terms  originated  by 
Southern  and  Goodinqs^  to  preserve  rotational  invariance  when  the  nonuniform 
strains  associated  with  finite-Wr-  /elength  elastic  waves  are  considered.  The 
additional  consideration  of  the  quadratic  nature  of  finite  strains,  first 
considered  in  the  context  of  dynamic  magnetoelastic  effects  in  rare-earth  materials 
by  Southern  and  Goodings,  is  included  as  well,  but  many  of  the  effects  treated  here 
are  not  strongly  dependent  on  the  nonlinearity  of  the  lattice  strain. 

The  principal  contributions  to  the  model  Hamiltonian  are  summarized  in  the 
following  paragraphs: 

1 .  Exchange  Interaction:  Molecular-Field  Approximation 

For  the  purposes  of  making  relatively  straightforward  calculations, 

the  usual  molecular-field  approximation  is  used  here  ~  represent  the  exchange 

interaction.  Although  this  approximatioi  is  subject  to  its  usual  limitations 

and  is  probably  not  a  very  good  approximation  at  very  low  temperatures,  it  is 

probably  quite  satisfactory  for  paramagnetic  materials,  which  are  the  subject  of 

some  of  the  experimental  work  described  in  subsequent  sections  of  this  report, 

and  the  molecular-field  approximation  seems  to  give  satisfactory  results  in  several 

cases  in  which  rare-earth  materials  are  considered  near  their  magnetic  phase 

transitions.  The  fact  that  the  pure  rare-earth  elements  and  many  of  their  alloys 

and  compounds  exhibit  periodic  ordering  structures  can  still  be  taken  into  account 

with  the  molecular-field  approximation,  if,  as  shown  by  many  authors,  particu- 
13 

larly  Nagamiya  ,  the  effective  molecular  field  is  taken  to  be  spatially  periodic 
in  accordance  with  the  periodic  nature  of  '.he  magnetization.  Although  it  shown 
by  several  authors^ ' that  there  are  sigi  ificant  anisotropic  contributions  to 
the  exchange  interaction  in  rare-earth  materials,  a  prediction  verified  by 
experiments  on  inelastic  neutron  scattering^ ,  the  treatment  given  here  does 
not  take  the  anisotropic  part  of  the  exchange  interaction  into  account  except  ’n 
the  case  of  magnetoelastic  interactions,  as  discussed  below.  This  approximation 
does  not,  apparently,  impair  greatly  the  validity  of  most  of  the  results  described 
here . 

That  part  of  the  model  Hamiltonian  representing  the  static  exchange  interac¬ 
tion  can,  therefore,  be  represented  in  the  following  way: 
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In  this  expression,  g  is  taken  to  be  the  free-ion  Lande  factor,  is  the  Pohr 

-  ►  0  |  ( 

maqneton,  H  is  the  effective  molecular  field  acting  uniformly  on  all  atoms  of 

th  ^ 

thn  ri  plane  of  ions  perpendicular  to  the  periodic  orderim  wave  vector,  and 

J.  is  the  angular-momentum  operator  for  the  i^  ion  of  the  n^  plane.  If 
m 

there  is  no  periodicity  of  finite  wavelength  in  the  magnetic  ordering,  then,  of 
course,  the  effective  molecular  field  is  uniform.  When  periodic  orderinq  is 
present,  however,  the  effective  molecular  field  will  exhibit  the  same  periodicity 
as  that  of  the  magnetization.  In  most  calculations,  the  molecular  field  will  be 
proportional  to  the  magnetization  or,  in  the  case  of  periodic  ordering,  a 
function  of  the  magnetic  moments  of  a  number  of  planes  of  ions.  Its  specific 
nature  must  be  considered  for  each  type  of  calculation  performed. 

In  most  of  what  follows,  better  approximations  than  the  molecular-field 
approximation  could  undoubtedly  be  employed,  although  with  a  great  increase  in 
the  difficulty  of  performing  detailed  calculations.  In  particular,  for  calcula¬ 
tions  involving  elastic  waves  of  very  high  frequencies  at  temperatures  near  the 
magnetic  ordering  transition  temperature,  the  random-phase  approximation  would 
provide  more  reliable  results,  whereas  for  low-temperature  calculations  a  fully 
quantized  treatment  of  the  coupled  spin-wave-phonon  system  would  be  necessary. 

For  most  of  the  problems  of  irterest  here,  however,  the  frequency  of  interest 
is  relatively  low  (10-700  MHz) ,  and  resonant  spin-wave-phonon  interactions 
do  not  seem  to  be  important  for  the  determination  of  the  characteristics  of  most 
of  the  experimentally  observed  phenomena  reported  here. 

2 .  Anisotropy  Energy:  Crystal-Field  Approximation 

Although  anisotropic  exchange  interactions  undoubtedly  constitute  a 

,  ■  _  14-16 

source  of  magnetic  anisotropy  energy  m  many  rare-earth  materials  ,  the 

static  anisotropy  will  be  t'eated,  following  the  example  of  many  previous  authors1, 
in  terms  of  an  effective  Hamiltonian  operator  representing  the  interaction  of  each 
individual  rare-earth  ion  with  the  electric  field  produced  by  its  environment. 

This  crystal-field  approximation  is  undoubtedly  quite  satisfactory,  since  the  f- 
electrons  of  the  rare-earth  ions  are  highly  localized,  but  a  full  explanation  of 
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certain  observed  magnetic  o^deiing  structures  and  certain  features  of  the  static 
magnetostriction  of  rare-earth  element'  may  require  the  introduction  of  the  aniso¬ 
tropic  exchange  interaction  omitted  here.  Most  of  the  dynamic  magnetoelastic  effects 
of  interest  here,  however,  appear  to  be  explained  adequately  without  the  inclusion 
of  anisotropic  exchange  in  the  static,  equilibrium,  terms  of  the  model  Hamiltonian. 

^he  model  Hamiltonian  representing  the  crystal-field  single-ion  anisotropy 
can  be  expressed  in  terms  of  angular-momentum  operator  equivalents,  through  the 
use  of  the  Wigner-Eckart  theorem,  as  first  shown  by  Stevens^  .  The  angular-momen¬ 
tum  tensor  operators  employed  here  are  essentially  the  same  as  those  used  by 
different  authors,  but  the  normalization  factors  have  been  chosen,  in  some  cases, 
differently  in  order  to  facilitate  numerical  calculation.  The  exact  definition 
of  the  operators  used  here  is  given  in  Appendix  I.  For  hexagonal  symmetry,  such 
as  that  of  the  hep  rare-earth  elements,  the  appropriate  model  Hamiltonian  is  the 
following : 


-I 


P°(Jon)  . 

2  20  l 


60*  i 


+  l 
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(II-2) 


In  this  exnression ,  the  coefficients  P^  depend  on  the  nature  of  the  crystalline 
environment  and,  in  principle,  can  be  calculated.  However,  since  most  reported 
calculated  values  are  not  in  particularly  good  agreement  with  experimentally 
determined  values  for  these  coefficients  for  most  rare-earth  materials  ,  onlv 
the  values  estimated  from  experimental  measurements  are  u  ied  here.  In  eq.  (II-2) 
the  summation  is  taken  over  all  ions  in  the  crystal,  and  it  is  assumed  that  each 
ion  interacts  with  the  lattice  in  the  same  way,  so  that  the  coefficients  are  identic 
for  all  ions. 

Since  the  anisotropy  free  energy  is  an  appreciable  fraction  of  the  exchange 
energy  for  many  rare-earth  materials,  it  must  be  taken  into  account  explicitly  in 
most  calculations  of  the  magnetic  and  r.agnetoelastic  properties.  For  many  purposes, 
however,  the  anisotropy  energy  can  be  represented  adequately  by  means  of  only  the 
first  term  of  eq.  (II-2) ,  the  second-degree  uniaxial  anisotropy  term,  and  this 
truncation  of  the  complete  anisotropy  model  Hamiltonian  will  be  carried  out  in  most 
of  the  work  described  here,  for  the  sake  of  the  resulting  simplification  of  the 
calculations . 
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3 •  MagneLoolastic  Interactions 

Two  types  of  magnetoelastic  interaction  have  been  considered  by 
a  number  of  authors,  following  the  general  formalism' of  Callen  and  Callen7. 

The  first  type,  which  is  designated  tne  single-ion  interaction,  results  from 
the  strain  dependence  of  the  crystal  field.  This  strain  dependence  produces 
as  one  effect  a  modification  of  the  coefficients  of  the  anisotropy  terms  of  the 
model  Hamiltonian  described  above,  but  since  certain  strain  compo„er-s  reduce  the 
point-group  symmetry  at  each  lattice  site,  they  also  introduce  crystal-field 
terms  of  lower  symmetry  than  those  characteristic  of  the  unstrained  lattice. 

The  second  type  of  magnetoelastic  interaction  is  the  so-called  two-ion  interaction, 
which  results  from  the  strain  dependence  of  the  exchange  interaction,  m  the 
usua.  description  of  this  interaction7  it  is  customary  to  include  the  effects  of 
anisotropic  exchange  interactions  although,  as  discussed  above,  they  are  often  not 
taken  explicitly  into  account  in  the  equilibrium  exchange  Hamiltonian  of  eq.  (11-1) , 
Both  of  these  types  of  magnetoelastic  interaction  are  considered  h-re, 
although  some  of  the  calculations  make  explicit  use  of  only  the  single-ion' 
magnetoelastic  interaction  for  the  sake  of  simplifying  calculations.  In  most 
oases  both  the  single-ion  and  two-ion  interactions  lead  to  qualitatively  similar 
observable  effects  such  that  it  would  be  difficult  to  separate  the  two  interac¬ 
tions  on  the  basis  of  experimental  measurements.  Following  Southern  and  Goodings11 
the  treatment  presented  here  includes  both  the  nonlinear  nature  of  finite  strains 
and  the  added  effects  on  the  behavior  of  shear  elastic  waves  resulting  from  the 
rotational  invariance  imposed  on  the  magneto,  lastic  interaction  by  these  authors. 

a-  single-Ion  Magnetoclastic  Intorachi nn 


There  are  two  principal  contributions  to  the  single-ion  magneto¬ 
elastic  interaction.  The  first  is  due  to  the  strain  dependence  of  the  crystal  electric 
field  resulting  from  the  lattice  deformation  due  to  the  strain.  This  contribution 
is  the  one  normally  treated  in  problems  involving  static  magnetostriction.  The 
second  contribution  arises  from  the  effective  rotation  of  the  crystal  field 
resulting  from  the  strain  associated  with  shear  elastic  waves  of  finite  wavelength 
It  is  this  second  contribution  which  was  introduced  by  Southern  and  Goodings11  as 
a  result  of  their  imposition  of  rotational  invariance  on  the  model  Hamiltonian, 
following  the  work  of  Meloher  and  Brown20.  Each  of  these  contributions  is  con- 
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sidered  separately  in  what  follows. 

The  normal  single-ion  magnetoelastic  interdiction  due  to  deformation  of  the 
crystal  lattice  as  the  result  of  a  homogeneous  strain  would  be  expecued  to  contain 
contributions  due  to  the  strain  dependence  of  crysta] -field  terms  of  degree 
1=2,  4,  and  6,  including  all  appropriate  m-values.  For  the  sake  of  simplifica¬ 
tion,  however,  only  the  terms  of  lowest  degree,  with  l  =  2,  are  included  here. 

It  would  also  be  expected  that,  for  each  value  of  £,,  terms  linear  in  the  lattice 
strain  and  higher-order  strain  terms  should  contribute  to  the  interaction.  In 
fact,  Southern  and  Goodings  took  account  of  a  part  o^  the  higher-order  strain 
terms  by  using  the  conventional  definition  of  finite  lattice  strain,  which 
includes  terms  of  second  degree  in  the  particle-displacement  gradient.  That  is, 
the  conventional  definition  of  strain  is  expressed  in  the  following  wa^ : 


.  3u.  3u.  .  3u,  3u. 

1_  .  i  1_  X  A 

ij  2  3x  3x.  +  2  3x.  3x. 

j  i  13 


(II-3) 


In  this  expression,  the  u^  represent  the  components  of  the  vector  displacement 
of  any  lattice  point  from  its  equilibrium,  position.  The  first  term,  in  brackets, 
represents  the  usual  small-signal  strain  tensor.  This  definition  of  strain  was 
originally  chosen  to  satisfy  the  requirement  that  the  strain  tensor  represent  a 
pure  deformation  of  the  lattice,  but  it  is  not,  of  course,  the  only  tensor  which  car. 
be  composed  of  linear  and  quadratic  combinations  of  the  displacement  gradient. 

In  fact,  such  a  combination  would  not  manifestly  arise  in  an  expansion  of  the 
crystal  field  in  terms  of  the  displacement  gradient  (nor  does  it  arise  naturally 
in  the  usual  expansion  of  the  crystal  lattice  energy  in  terms  of  second-  and  third- 
order  elastic  constants).  However,  the  strain,  as  defined  in  eq.  (II-3) ,  does 
provide  a  useful  symmetric  tensor  for  crystal-field  and  other  expansions,  with  the 
advantage  that  its  direct  product  with  itself  is  also  a  tensoi  of  higher  ran... 
Consequently,  the  form  of  eq.  (II-3)  is  used  here,  although  in  many  cases  only 
the  small-signal,  linear  part  of  the  strain  need  be  used  in  a  specific  calculation. 
Wren  only  the  small-signal  part  is  used,  it  will  be  designated  by  the  svmbol  e  . 

In  addition,  it  is  often  convenient  to  use  the  six-dimensional  Voigt  notation  for 
such  quantities  as  strain  and  the  elastic  constants.  For  the  strain  (small-signal 
values  only),  the  notation  e^  will  be  used  where  the  normal  convention  (e^  =  , 

e2  =  r22'  e3  =  G33 '  e4  =  G23  +  C32'  e5  =  C13  +  Si'  e6  =  C12  +  Sl}  wil1  be  followed- 
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Ihe  following  expression  will  be  used  as  the  model  Hamiltonian  for  that 
part  of  the  magnetoelastic  interaction  due  to  single-ion  interactions  with  the 
crystal  field,  not  including  that  part  due  to  the  rotation  of  the  crystal  field, 
which  is  treated  separately  below: 
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me 


c  b“'2  eJ  *  b-  e«i  i  u2oh  ♦  b^# 2 1  ,Ej ♦ 
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(II-4) 


In  this  expression,  which  is  valid  for  the  hexagonal  symmetry  wnifh  holds  for  the 
rare-earth  elements,  lines.-  combinations  of  the  cartesian  strains  which  form  the 
components  of  the  various  irreducible  representations  of  the  hexagonal-symmetry 
point  group  appropriate  to  the  hep  lattice  have  been  employed.  The  definitions  of 
these  quantities,  denoted  by  the  symbol  e“,  which  differ  slightly  from  those 
employed  by  Callen  and  Callen7,  are  given  m  Appendix  II.  It  should  be  remembered 
that  the  magnetoelastic  model  Hamiltonian  given  in  eq.  (II-4)  ignores,  for  the 
sake  of  simplification,  higher-order  terms  in  both  the  strain  and  the  a.igular- 
momentum  tensor  operators.  These  lowest-order  terms  are,  however,  sufficient  for 
the  description  of  most  magnetoelastic  phenomena  in  rare-earth  materials. 

That  part  of  the  magnetoelastic  interaction  due  to  the  rotation  of  the  crystal 
field  associated  with  shear  strain  waves  of  finite  wavelength  may  be  treated  to 
the  same  order  in  the  following  way.  If  a  static  external  stress  were  applied  to 
a  crystal,  of  course,  only  symmetric  strains,  both  shear  and  longitudinal,  would 
result,  A  pure-mode  shear  elastic  wave  propagating  through  a  crystal  with  finite 
wavelength,  however,  is  characterized  by  a  skew-symmetric  particle-displacement 
gradient,  and  the  resulting  motion  of  the  lattice  may  be  described  in  terns  of 
both  a  symmetric  (small-signal)  strain  wave  and  an  antisymmetric  rotation  of  the 
coordinate  system.  For  example,  a  shear  wave  propagating  along  the  z-axis  with  its 
polarization  along  the  x-ax„  would  be  characterized  by  a  particle-displacement 
gradient  Su^/Sz,  which  coaid  be  written  as 


&u 


(II-4) 


Thus,  following  Southern  and  Goodings^,  a  general  elastic  wave  cculd  be 
represented  by  a  generalized  strain,  including  both  symmetric  and  antisymmetric 
parts,  according  to  the  following  notation: 
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The  symmetric  part  of  the  generalized  expression  in  eq.  (II-6)  is  just  the 
small-signal  strain,  which  is,  of  course,  treated  in  the  expression  of  eq.  (II-4) 
The  antisymmetric  part  is,  however,  equivalent  to  a  rotation  of  the  coordinate 
axes,  and  this  rotation  leads  to  a  term  in  the  model  Hamiltonian  for  the  single¬ 
ion  magnetoelastic  interaction  which  reflects  the  effect  of  the  rotation  on  the 
anisotropy  terms  given  in  eq.  (II-2)  .  If  only  the  second-degree  anisotropy  term 
is  considered,  for  the  sake  of  simplification  and  because  the  fourth-  and 
sixth-degree  terms  are  not  included  in  eq.  (II-4)  ,  then  the  following  terms 
mus-  be  added  to  the  single-ion  magnetoelastic  interaction: 


Hme  2  P2  ^  tfi13(J21)i  +^\3(J21)i]  *  (H-7) 

In  the  above  expression,  the  quantities  and  are  defined  as  follows: 
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As  pointed  out  by  Southern  and  Goodings,  this  contribution  to  the  single-ion 
magnetoelastic  interaction  should  lead  to  observable  effects  in  the  consideration 
of  the  magnetoelastic  contributions  to  those  elastic  constants  associated  with 
shear  strain  components. 


b .  Two-Ion  (Exchange)  Magnetoelastic  Interaction 


The  two-ion  magnetoelastic  interaction  arises  from  the  strain 
dependence  of  the  exchange  interaction,  and  it  can  be  written  down  to  first 
order  in  the  strain  in  tho  following  way; 


(II-9) 
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It  should  be  noted  that  the  expression  for  the  two-ion  magnetoelastic  interaction 
given  in  eq.  (11-10)  differs  from  those  used  by  Callen  and  Callen7  and  by  Southern 
and  Goodings  because  of  a  different  normalization  of  the  strain  and  angular- 
momentum  operators  employed  here.  It  might  be  argued  that  an  additional  term 
due  to  the  rotation  associated  with  finite-wavelength  shear  strain  waves  should 
also  be  considered  here,  as  in  the  case  of  the  single-ion  magnetoelastic  interactio 
However,  the  dominant  contribution  to  the  static  exchange  interaction  in  the  unstra 
lattice  is  the  isotropic  Heisenberg  interaction,  and  any  anisotropic  equilibrium 
exchange  interactions  have  been  neglected  in  this  treatment.  Only  such  anisotropic 
contributions  to  the  exchange  interaction  would  lead  to  rotational  effects  in  the 
presence  of  shear  strain  waves. 

The  two-ion  magnetoelastic  interaction  is  easily  treated  within  the  frame¬ 
work  of  the  molecular-field  model,  and  it  is  included  where  necessary  in  the 
detailed  calculations  presented  below.  However,  this  interaction  usually  does  not 
lead  to  any  significant  qualitative  difference  from  results  obtained  through  the 
use  of  the  single-ion  interaction  alone,  and,  consequently,  it  need  not  always  be 
included  explicitly  in  specific  calculations.  Its  importance  in  each  case 
will  be  considered  in  what  follows. 

4.  Elastic  Energy — Long-Wavelength  Classical  Limit 

The  problems  of  interest  here  concern  the  generation  and  propagation  of 
coherent  elastic  waves  of  frequencies  from  10-20  MHz  up  to  the  X-band  microwave 
range  (9-10  GHz).  Over  this  entire  frequency  range,  the  wavelength  is  long  enough 
that  the  crystal  will  appear  to  be  a  continuum  insofar  as  its  effects  on  the 
wave  velocity  and  attenuation  are  concerned.  Thus,  the  purely  elastic  part  of 
the  problem  is  treated  here  in  the  low-frequency  continuum  limit,  while  the 
magnetic  and  magnetoelastic  parts  of  the  problem  are  treated  quantum-mechanically , 
in  the  same  sense  as  the  usual  semiclassical  treatment  of  the  interaction  of 
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radiation  with  natter.  In  this  approximation,  the  elastic  potential  energy  may 
be  written  in  the  usual  manner: 


u  =  7  I 


6  2  ijkl  Ci^lEijEkl  +  3  . 


—  y 

ijklmn  i3klmn  ij  bl 


E  +  •  •  • 

mn 


(11-11) 


In  this  expression,  the  elastic  energy  is  written  in  terms  of  the  elastic 
finite  strains  discussed  above.  An  exactly  equivalent  expression  may  be  written 
in  terms  of  the  small-signal  strains,  with  the  appropriate  c.  .  replacing  the 
finite-strain  E...  The  only  difference  between  the  two  ways  of  writing  the 
elastic  energy  is  the  fact  that  the  third-order  elastic  constants  c.  .  will 
be  different  in  the  two  representations  (as  would,  of  course,  any  elastic 
constants  of  still  higher  order) .  it  should  be  noted  that  the  second-order 
elastic  constants  may  be  defined  from  the  above  expression  for  the  elastic  energy 
as  the  second  partial  derivatives  of  the  energy  with  respect  to  the  strain  com¬ 
ponents  (in  either  the  finite-strain  or  small-signal  representations: 


32U 


32U 


°ijkl  3E .  .  3E 

i]  kl 


3e. . 3e  , 
13  kl 


(IJ-12) 


This  representation  for  the  elastic  constants  as  the  second  partial  derivatives 

elastic  energy  with  respect  to  the  strain-tensor  components  leads  natural¬ 
ly  to  the  possibility  of  extending  this  result  to  include  not  only  the  elastic 
energy  itself,  but  also  the  internal  energy  associated  with  the  magnetic  and 
magnetoelastic  Hamiltonians.  For  the  case  of  low-frequency  elastic  strain,  . 
this  approach  would  lead  to  a  simple  method  for  the  calculation  of  the  magneto¬ 
elastic  contributions  to  the  elastic  constants.  Since  such  magnetoelastic  contri¬ 
butions  to  the  elastic  constants  of  rare-earth  materials  are  one  of  the  most 
readily  observable  of  the  dynamic  magnetoelastic  phenomena,  as  seen  in  the  followino 
sections  of  this  report,  this  point  will  be  explored  at  some  length  below. 

The  kinetic  elastic  energy,  in  the  classicil  limit,  is  given  as  the  volume 

integral  of  the  classical  kinetic  energy  of  each  volume  element  of  the  lattice 
in  the  usual  manner 
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B.  Method  of  Calculation 


Because  each  type  of  rare-earth  material  presents  a  different  physical 
configuration  of  the  ions,  leading  to  different  detailed  computational  procedures 
for  each  of  the  many  possible  materials,  detailed  results  of  calculations  based 
on  the  theoretical  model  described  above  are  not  give  here.  The  general  method 
of  approach  is  described  briefly,  and  specific  results  are  given  in  the  following 
sections,  m  which  the  experimental  program  and  its  results  are  described. 

If  the  model  described  in  the  preceding  section  is  adopted,  the  state  of  a 
crystalline  specimen  as  a  whole  is  described  in  terms  of  the  thermodynamic  functions 
which  can  be  calculated  on  the  basis  of  the  model  Hamiltonian.  For  example,  in 
the  calculation  of  the  magnetoelastic  contributions  to  the  elastic  constants, 
it  is  necessary  to  compute  the  internal  energy,  whose  second  partial  derivatives 
are  the  conventionally  defined  adiabatic  elastic  constants.  The  calculation  of 
a  state  function  such  as  the  internal  energy,  even  in  a  simplified  model  such 
as  that  described  here,  requires  rather  complicated  numerical  computation  even 
for  paramagnetic  or  simple  ferromagnetic  materials.  When  periodic  ordering  is 
also  considered,  the  computations  become  even  more  complicated.  As  a  result,  only 
a  very  few  special  cases  have  been  worked  out  in  any  detail. 

The  general  approach  is  the  calculation  of  a  quantum-mechanical  density 
matrix  operator.  Then,  the  state  variables  may  all  be  calculated  by  means  of 
evaluating  straightforward  partition  sums.  Any  partial  derivatives  which  are  then 
needed  for  the  determination  of  elastic,  magnetic,  or  magnetoelastic  properties 
may  usually  be  worked  out  numerically  or,  in  some  cases,  analytically.  The 
problem  is  complicated,  however,  in  the  case  of  the  elastic  constants  and  other 
second-order  properties,  requiring  the  calculation  of  second  derivatives  of 
state  variables,  in  that  the  partition  function  or  density  matrix  may  vary  in 
first  order  with  respect  to  the  strain  of  other  variable  of  interest.  Thus,  the 
actual  computation  of  magnetic  and  magnetoelastic  properties  may  be  rather 
difficult  in  many  cases.  Those  specific  results  which  have  been  worked  out 

within  the  framework  of  this  theoretical  approach  are  described  in  the  following 
sections. 

It  should  be  pointed  out  that  the  approach  taken  here  is  probably  appropriate 
only  when  elastic  waves  whose  wavelength  is  long  compared  to  lattice  spacings,  and 
whose  frequencies  are  low  compared  to  natural  relaxation  rates,  are  considered. 
Otherwise,  r  coupled  spinwave-phonon  approach  using  second-quantization  techniques 


m 
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would  be  required.  Such  an  approach  was  employed  by  Southern  and  Goodings11, 
but,  for  the  reasons  mentioned  above,  and  because  of  its  greater  complexity  and 
difficulty  of  computation,  this  method  was  not  employed  in  the  theoretical 
interpretation  of  the  experimental  work  described  below. 

III.  EXPERIMENTAL  OBSERVATIONS  OF  DYNAMIC  MAGNETOELASTIC  EFFECTS  IN  RARE-EARTH 
MATERIALS  AND  THEIR  INTERPRETATION. 

This  section  gives  first,  an  overview  of  dynamic  magnetoelastic  effects  in 
rare-earth  materials,  Then  a  description  of  work  carried  out  in  this  program  on 
magnetoelastic  effects  in  terbium,  erbium,  and  in  various  rare-eaith-iron-group 
materials  is  given. 


A*  £ynamic  Magnetoelastic  Effects  in  Rare-Earth  Materials:  A  Survey 


The  extremely  strong  magnetoelastic  interactions  found  in  many  rare-earth 
materials  lead  to  a  number  of  interesting  observable  phenomena,  the  most  striking 
of  which  is  the  so-called  "giant"  magnetostriction  exhibited  at  low  temperatures 
by  the  single  crystals  of  the  pure  elements  Tb,  Dy,  Ho,  and  Er18.  Recently,  C" ark 
and  Be Ison  have  reported  that  a  new  class  of  compounds,  the  cubic  Laves-phas* 
materials,  RFe2,  where  R  represents  any  of  several  rare-earth  elements,  also 
exhibit  very  large  magnetostriction,  even  in  polycrystalline  specimens  at  room 
temperature.  The  static-equilibrium  magnetostriction  of  the  pure  rare-earth  elements 
is  at  present  fairly  well  understood  in  terms  of  a  theory  due  to  Callen  and  Callen7, 
who  took  into  account  both  single-ion  magnetoelastic  interactions  arising  from  the 
strain  dependence  of  the  the  crystal-field  anisotropy  energy  and  two-ion  magneto¬ 
elastic  interactions  arising  from  the  strain  dependence  of  the  exchange  interaction 
energy  between  pairs  of  ions.  The  experimental  results  of  Clark,  DeSavage,  and 
Bozorth  on  the  static  magnetostriction  of  single-crystal  dysprosium  are  in  parti¬ 


cularly  good  agreement  with  the  theory  of  Callen  and  Callen. 

This  report  is  concerned  with  magnetoelastic  interactions  in  situations  invol¬ 
ving  time-dependent  elastic  strains.  Two  principal  areas  are  of  interest  here.  The 
first  is  the  effect  of  magnetoelastic  interactions  on  the  propagation  of  elastic 
waves  in  rare-earth  materials.  Both  the  velocity  and  the  attenuation  of  elastic 
waves  in  magnetoelastic  materials  exhibit  anomalous  behavior  compared  to  nonmagnetic 
materials,  and  the  effect  of  the  magnetoelastic  interactions  on  the  elastic-wave 
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velocity  has  been  Interpreted  primarily  in  terns  of  the  adiabatic  elastic  constants 
of  rare-earth  materials.  These  elastic  constants  are  usually  measured  by  measuring 
the  velocities  of  elastic  waves  in  the  frequency  range  of  several  Mils  along  certain 
crystallographic  directions,  and,  in  the  case  of  the- rare-earth  elements,  the 
temperature  dependences  of  the  elastic  constants  exhibit  behavior  which  is  anomalous 
compared  to  that  of  most  other  materials.  The  temperature  dependence  of  the  elastic 
constants  of  several  rare-earth  elements  is  described  in  what  follows,  and  an 
explanation  of  the  magnetoelastic  contributions  to  the  elastic  constants  is  given 

in  terns  of  an  extension  of  the  theory  of  Callen  and  Callen  to  the  domain  of  time- 
dependent  lattice  strain. 

The  second  area  of  interest  is  the  magnetoelastic  generation  of  elastic  waves. 
The  phenomenon  of  magnetostriction  has  been  used  for  many  years  in  ultrasonic 
transducers,  with^many  technological  applications,  when  it  was  found,  by  Legvold, 
Alstcd,  and  Rhyne  that  single- crystal  dysprosium  and  holmium  both  exhibit  magneto¬ 
striction  at  low  temperatures  some  100  times  greater  than  that  of  nickel  and  other 
iron-group  materials,  the  interesting  possibility  was  raised  that  such  materials 
might  lead  to  the  development  of  highly  efficient  ultrasonic  transducers,  perhaps 
capable  of  operation  at  very  high  frequencies.  Haley,  Donoho,  and  Blackstoad24 
reported  ultrasonic  elastic-wave  generation  in  thin-film  polycrystalline  trans¬ 
ducers  of  Gd,  Dy,  Ho,  and  Er  at  frequencies  between  1,000  and  10,000  MHz,  with 
efficiency  somewhat  greater  than  that  observed  in  similar  thin-film  „i  transducers 
The  efficiency  of  the  rare-earth  transducers  was  not,  however,  as  high  as  would 
have  been  expected  from  a  simple  comparison  of  the  static  magnetostriction  of 
these  rare-earth  elements  with  that  of  nickel.  This  deficiency  was  attributed  in 
part  to  the  lack  of  purity  of  the  films  and  to  their  polycrystalline  nature. 

The  rare-earth  transducers  also  possessed  the  serious  defect,  from  the  standpoint 

of  using  them  in  practical  applications,  that  they  required  low  temperatures  and 
high  magnetic  fields  for  their  operation. 

The  recent  discovery  by  Clark  and  Belson21  that  the  intermetallic  compound 
TbFe2  exhibits  very  large  static  magnetostriction  at  room  temperature  in  polycrys- 
tallinc  specimens  suggested  that  this  new  material  might  avoid  the  difficulties 
mentioned  above  and  lead  to  the  possibility  of  the  development  of  worthwhile  new 
magnetostrictive  transducers.  Magnetostrictive  ultrasonic  generation  was  observed 
in  thin  films  of  this  material  by  Donoho  et  al25  at  room  temperature  and  at  fre¬ 
quencies  from  300  MHz  to  1,500  MHz.  The  results  which  have  been  obtained  to  date 
on  magnetostrictive  generation  in  the  pure  rare-earth  elements  and  in  Tb-Fo  films 
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arc  reviewed  briefly  in  this  section,  and  a  more  detailed  account  of  the  results 
obtained  in  this  program  on  Tb-Fe  compounds  is  given  in  a  subsequent  section. 

1*  M-a9net°elastic  Contributions  to  Elastic  Constants 

The  temperature  dependences  of  certain  of  the  elastic  constants  of  Tb,  Ho, 

and  Er  are  given  in  Figures  2-  10.  The  elastic  constants  of  Tb  and  Ho  were 
measured  by  Salama  otal  '  27,  and  those  of  Er  have  heen  measured  by  HubbeU 

—  “29‘  In  addition'  elastic  constants  of  Gd  have  been  measured  by  Long 
et  al  ,  and  those  of  Dy  have  been  measured  by  Rosen  and  Klimker30.  However, 
only  results  obtained  in  the  present  program  are  shown  in  the  figures.  For  I 
normal,  nonmagnetoelastic,  material,  the  elastic  constants  generally  increase 
smoothly  with  decreasing  temperature  as  the  material  becomes  elastically  "harder-. 
As  seen  in  Figs.  2-10,  however,  there  are  numerous  departures  from  this  normal 
behavior  in  the  case  of  the  rare-earth  elements. 

in  the  case  of  Tb  as  the  temperature  is  lowered,  there  is  a  paramagnetic-helimag- 

P  .o  transition  at  _  ..30  K  and  a  he]  magnetic- ferromagnetic  transition 

at  Tc  -  219  K.  At  both  of  these  transition  temperatures  the  elastic  constants  of 

Tb  exhibit  anomalous  behavior.  The  elastic  mnet-anh  _  f 

elastic  constant  c^,  for  example,  has  cusp- 

like  anomalies  at  these  temperatures,  whereas  cJ3  exhibits  even  more  pronounced 

anomalies,  as  does  c^.  On  the  other  hand,  c„  undergoes  only  a  change  in  slope  a, 

the  second-order  transition  at  V  while  it  exhibits  a  small  cusp  at  the  first 
order  transition  at  T  . 


c„? 
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Fig.  2  Temperature  dependence  of  cn  for  single-crystal  Tb. 
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Fig.  3  Temperature  dependence  of  c33  for  single-crystal 


Tb. 


Fig.  4  Temperature  dependences  of  c^  and  c&6  for  single-crystal  Tb. 

For  the  elastic  instants  of  Ho,  the  qualitative  nature  of  the  magnetoelas¬ 
tic  contributions  to  the  elastic  constants  is  similar,  as  seen  in  Figs.  4-7.  in  this 
case,  unfortunately,  strong  elastic-wave  attenuation  prevented  measurements  below 
approximately  80  K,  and  only  the  anomalies  at  the  second-order  paramagnetic- 
helimagnetic  phase  transition  at  T^  =  130  K  can  be  seen.  The  first-order  helimag- 
netic-ferromagnetic  transition  at  Tc  =  20  K  could  not  be  observed.  In  the 


case  of  Ho, 
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it  can  be  seen  that  e,,  undergoes  only  a  change  in  slope  at  T  ,  but  that  c  exhibits 

chi:;i  ;:;1:r;:i;anomaiy  at  tMs  tmperature-  *>“  -  =e6  ~ 

h  ““  °f  erbl“"  is'  in  so™  th<=  «>»t  interesting  of  those  discussed 

Er  as  v,15  W°rt  lla  t0  C°nSlder  in  som!!  detail  the  magnetic  ordering  phases  of 
Er,  as  shown  m  Fig.  n# 


82-  'olmium 


Temperature  [  k] 


Fig.  5  Temperance  dependence  of  c  for  single- 


-crystal  Ho. 
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Fig,  6  Temperature  dependence  of 


c33  for  single-crystal  Ho. 
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Fig. 


Temperature  dependence  of  c  for  , 

44  ror  single-crystal  Ho. 
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Pi9.  8  Temperature  dependence  .f  c66  for  single-crystal  Ho. 

T  .  “a  t  “t"'  ^  "  Para"a9netiC-  “h— =  «“  temperatures  between 
0  8  1  anti ferromagnetic,  with  a  sinusoidal  modulation  of  th 

component  of  the  magnetisation  of  each  basal  elan  ,  «  C’a> 

and  Tp  the  ordering  is  such  that  the  c-axis  oompon  Tf  l  10"S'  Bet“een  ^  =  2<> 
plane  of  ions  exhibits  a  square-wave-like  .  °  "  magnetization  of  each 

another,  while  the  basal  ol  Savior  in  going  £rom  one  plane  to 

basal-plane  component  is  helically  ordered.  Below  Tp  the  order! 
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fe  with  .  uniform  c-axis  component  of  the  magnetization 

e  ica  or  errng  of  the  basal-plane  component  of  the  magnetization. 


ig.  9  Temperature  dependences  of  c  and  c 

11  33 


for  single-crystal  Er. 


& 


Fig.  10  Temperature  dependences 


of  C44  and  c56  for  single-crystal  Er. 
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Fig.  11  Magnetic  Ordering  Phases  of  Single-Crystal  Er. 


irst-order  transition.  The  neutron-diffraction  data  of  Cable  et  a,31  do'not 
however,  permit  any  dedurU™  ^  —  —  do  not' 


- - uiuii  udta  or  Cable  et  a  1 

however  permit  any  deduction  to  be  made  concerning  the  transition^  at"  T  If  h 
transition  were  made  abrunt.lv  at  TQ  If  thls 


transit-inn  ,,  ,  —  —  rne  transition  at  T  If 

transition  were  made  abruptly  at-  t  Q 

this  temperature  as  tn„  ’  .  «  ‘°  structure  L 


this  temperature  as  the  temporal  was  lowered  from  a  ““ 

then  the  transition  would  surely  be  a  first  d  slnusordal  range  above  T 

if  the  helically  ordered  component  b  "  °"  “““  ^d, 

gradually,  starting  from  a  sere  ^^VlPT^01'  ^ 
zation  of  the  c-axis  magnetisation  at  this  ^  ^ 

be  of  second  order.  The  behavior  of  the  elast  •  ’  ransrtion  shoulc 

on  the  nature  of  thi-  phas  t  •  constants  at  T  sheds  some  loght 

Ptovided  by  the  be  a;  “on  -  also 

in  the  cases  of  all  the  LI  I  r  T  atten“ati°"  ^  ^  ^  "  “* 

viously  been  measured,  inclL  D y“  f “  “h°Se 

the  elastic  constant  ■  '  ^  “*  reS“ltS  «"  **  shown  here, 

transition  and  a  ousp-JL"  iLuLllT I  ^  “  ‘  ~ 

observation,  one  would  surmise  that  the  transition  at  tVTLT  1“  T 
transition.  The  evidence  from  tu  i  •  Q  '  ct'  a  second-ord< 

supports  this  conclusion  Lc  h  h  ^  ^  12  X: 

T  or  T  but  it  doe  shear-wave  attenuation  exhibits  no  peaks  at 

"  S  “  er,°  S  Steri>  ch“»*  at  the  first-order  transition  at  T 

C* 


ATTENUATION  (db/cm) 
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Fig. 


*-  Temperature  dependence 


of  longitudinal  c-axis  attenuation 


in  Er. 


ce  of  shear  c-axis 


Fig.  13  Temperature  dependen 


attenuation  in  Er. 
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.  -  -ir  “  r,“  “::r  r  ™ :  - — 

as  discussed  in  a  previous  section  of  this-  report  J.  h  7 

of  the  thermodynamic  internal  energy  of  the  aystll it ^  ^  dRriVatiVeS 
Of  the  lattice  ufrain  r  ^  respect  to  the  components 

-  -stic the  *> 

internal  energy  to  second  degree  in  Ule  straTflt  r0nPUtI!  T"  °f  the 

“r which 

arising  in  attempting  to  ■'  T 

success  nas  teen  obtained,  However,  in  describin,  H  b  b  ^ 

contributions  to  tbe  elastic  constants  of  ^ 

field  approximation  to  treat  the  exch*  •  .  '  ll21ng  a  m°l ocular- 

followed  by  Hubbell  et  al28  t  •  31196  lntCraCtl°nS '  A  similar  approach  was 

y  uooeii  et  al  to  explain  the  behavior  of  c  in  Fr  ,-r  «.u  ... 

of  the  Paramaynetic-antiforromagnetic  phase  transition  a”T  -  86  K  Th  V1ClnitY 

took  account  only  of  the  sinal^  •  N  86  K>  These  authors 

Previous  section  and  h  I  r'agnGt°elaStiC  ^^ion  defined  in  the 

oection,  and  they  used  results  obtained  by  Naganiya13  for  t-h 

Of  the  sinusoidally  nodulated  e-atis  magnetisation  below  T  in  t,  •  ^ 

the  magnetic  and  magnetoelastic  contributions  to  the  internal  en “  °n  °f 

Hubbell  et  al  were  able  internal  energy.  In  this  way, 

near  T  T~h  °bt!,ln  *  «t  to  the  behavior  of  c 

near  T^,  as  shown  in  Fig.  14.  33 

The  detailed  computation  of  the  temperature  dependence  of  all  the  , 
donstants  remains  a  difficult  urohle™  ■  •  11  *  elastlc 

magnetic-ordering  phases  of  these  materUlH  “T °f 

a  theory  of  the  magnetoelastic  TT.ZT'"  ^ 

earths  making  use  of  a  rotations  n  ■  ■  1C  constarits  of  the  rare 

treatments  of  dynamic  magnetoelastic'iZrMtiZ”^''- "et  lnClUded  “  PreVi°“S 
account  of  the  quadratic  nature  of  the  convention^  deli 72^'^  T  ^ 
strain,  was  developed  primarily  to  explain  the  dependence  of  the  .i“? 
on  applied  magnetic  fields,  and  it  should  permit  a  more  satisfarto  C°nSta"tS 

the  problem  of  the  temperature  dependence  of  the  elastic  constan  A„  ^ 

motivated  by  the  work  of  Southern  and  Goodings13  on  the  field  d  ' 

elastic  constants  of  paramagnetic  Tb  is  des  •„  a  •  dapendence  of  the 

gnetic  Tb  is  described  in  a  subsequent  section  of  this 


report. 
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Fig.  14  Comparison  of  experimentally  determined  temperature  dependence 

°f  c33  m  Er  near  Tg  =  86  K  with  theory  using  molecular-field  approximation  and 
single— ion  magnetoelastic  interaction. 

2-  Magnetos trictive  Generation  of  Ultrasonic  Waves 


The  magnetostrictive  generation  of  high-frequency  elastic  waves  in  rare- 
earth  materials  was  first  reported  by  Haley,  Blackstead,  and  Donoho24,  who  observed 
strong  magnetostrictive  elastic-wave  generation  at  frequencies  in  the  range 
1,000  to  10,000  MHz  in  films  of  2-ym  thickness  of  Gd,  Dy,  Ho,  and  Er  vacuum,  eva¬ 
porated  onto^substrate  delay  lines  of  quartz  and  sapphire.  This  work  was  extended  by 
Donoho  et  al  to  include  Tb  and  the  intemetallic  compound  TbFe2<  The  films  were 
placed  in  an  rf  cavity  in  a  region  of  maximum  rf  magnetic  field,  with  the  rf  field 
parallel  to  the  film  surface.  A  dc  magnetic  field  up  to  60  kOe  could  be  applied 
at  any  direction  with  respect  to  the  surface  of  the  film.  An  rf  pulse  was  applied  to 
excite  the  magnetostrictive  generation  of  the  elastic  waves,  which  then  gave  rise 
to  a  sequence  of  ultrasonic  pulse  echoes  in  the  delay-line  substrate.  Typical  results 
for  the  field  dependence  of  the  amplitude  of  elastic  waves  magnetostrictively  genera¬ 
ted  in  Dy  and  Tb  are  shown  in  Figs.  15  and  16. 

Two  important  features  of  the  magnetostrictive  generation  of  ultrasonic  waves 
m  the  pure  rare  ecrths  are  evident  in  Figs.  15  and  16.  First,  appreciable 
efficiency  of  generation  occurs  only  in  the  presence  of  large  applied  fields; 
second,  because  the  magnetic  ordering  temperatures  of  these  materials  are  in  the 
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cryogenic  range,  they  can  be  used  as  transducers  only  at  temperatures  so  low  that 
their  usefulness  in  practical  applications  is  severely  limited.  Another  somewhat 
disappointing  feature  of  the  magnetostrict.ive  generation  of  elastic  waves  in  pure 
rare-ea-th  films  is  the  fact  that  the  maximum  generation  efficiency  is  not  as 
great  as  might  have  been  expected  from  the  static  magnetostriction.  In  this  respect, 
the  polycrystalline  nature  of  the  films  used  to  study  ultrasonic  generation 
is  probably  responsible  for  the  low  efficiency. 


Fig.  15  Ultrasonic  echo  amplitude  for  Dy  film  on  X-cut  quartz  substrate  as 
a  function  of  applied  magnetic  field  perpendicular  to  plane  of  field. 


Fig.  16  Ultrasonic  echo  amplitude  for  Tb  film  on  X-cut  quartz  substrate  as 
function  of  applied  field 
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If  is  evident  that,  despite  the  enormous  static  magnetostriction  exhibited 
by  single  crystals  of  the  pure  rare-earth  elements  at  very  low  temperatures, 
these  materials  will  not  become  useful  as  magnetostrictive  transducers.  Conse¬ 
quently,  when  Clark  and  Belson21  reported  their  observation  that  the  intermetallic 
compound  TbFe.,  exhibits  very  large  magnetostriction  at  room  temperature  in  poly¬ 
crystalline  specimens,  it  appeared  that  this  new  material  would  probably  be  a  much 
better  candidate  for  the  development  of  a  good  high-efficiency  room-temperature 
ultrasonic  transducer.  Thin-film  specimens  of  TbFe2  were,  therefore,  deposited  on 
delay-line  substrates,  and  it  was  quickly  demonstrated25  that  ultrasonic  generation 
at  room  temperature  and  at  high  frequencies  (700  KHz)  was  feasible  with  this  material 
A  particularly  interesting  feature  of  the  magnetostrictive  ultrasonic  generation 
m  TbFe.,,  however,  and  a  totally  unexpected  feature,  was  the  fact  that  thin  films 
of  this  materials,  unlike  bulk  specimens,  exhibit  a  very  strong  coercive  force  and 
a  very  high  degree  of  anisotropy,  leading  to  a  strong  remanent  magnetization 
perpendicular  to  the  plane  of  the  film.  This  strong  remanence  permits  the  generation 
of  ultrasonic  waves  at  high  efficiency  without  the  requirement  of  an  applied  magnetic 
field.  The  results  obtained  on  Tb-Fe  films  of  various  compositions  are  given  in 
more  detail  in  a  subsequent  section  of  this- report. 

In  conclusion,  it  has  not  yet  been  successfully  demonstrated  that  rare-earth 
materials  will  compete  with  currently  available  ultrasonic  transducer  materials, 
but  the  compound  TbFe.,  shows  enough  promise  to  warrant  further  study. 

B*  Magnetoelastic  Coupling  in  Paramagnetic  Terbium 


As  discussed  above,  the  strong  magnetoelastic  interactions  found  in  most 
rare-earth  materials  lead  to  significant  magnetoelastic  contributions  to  the  elastic 
constants  these  materials.  The  principal  way  in  which  these  magnetoelastic 
contributions  are  manifested  is  in  the  anomalous  temperature  dependence  of  the 
elastic  constants,  particularly  near  magnetic  phase  transitions.  Such  anomalies 
have  been  observed  in  the  cases  of  Gd29,  Dy30,  Tb26,  „o27,  and  Er28,  and,  in  some 
cases,  they  amount  to  anomalous  deviations  of  the  elastic  constants  of  several 
per  cent  from  what  would  be  expected  for  a  non-magnetoelastic  material.  Attempts 
to  explain  these  magnetoelastic  contributions  to  the  elastic  constants  in  terms  of 
the  fundamental  magnetoelastic  interactions  of  the  rare-earth  elements  have,  however 
been  only  partially  successful  because  of  the  difficulty  encountered  in  calcula¬ 
ting  such  thermodynamic  functions  as  the  internal  energy  near  a  phase  transition 
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in  a  magnetically  ordered  material.  Southern  and  Goodings11  have,  however,  in  a 
paper  oonoerned  with  the  general  problem  of  the  magnetoelastic  contributions  to 
the  eiastic  properties  of  rare-earth  materials,  pointed  out  that  these  magnetoelastic 
n  u  roes  to  the  elastic  constants  would,  perhaps,  be  significant  even  in  the 
paramagnetic  temperature  range,  where  calculations  may  usually  be  made  with  greater 
confidence  than  in  the  magnetically  ordered  phases.  In  their  approach  to  this  prob- 
em,  southern  and  Goodings,  basing  their  work  in  part  on  earlier  work  of  Melcher19 
and  Brown  ,  developed  a  rotationally  invariant  theory  of  magnetoelastic  interactions, 
and  they  took  account  of  the  nonlinear  nature  of  finite  strains  to  arrive  at  a 
simple  result  for  rare-earth  materials  in  the  paramagnetic  phase. 

Motivated  by  the  work  of  Southern  and  Goodings11  in  calculating  the  magneto¬ 
elastic  contributions  to  the  elastic  constants  c  and  c,_  for  rare-earth 
mate  1S  in  the  paramagnetic  phase,  „  experi„e„t  _  ^  ^  ^ 

o  these  two  elastic  constants  on  temperature  and  on  the  value  of  an  applied 
magnetic  field  for  the  case  of  paramagnetic  terbium  at  temperatures  from  just 
above  Tn  -  229  K  to  approximately  MO  K.  The  results  of  this  experimental  work 
are  in  at  least  partial  agreement  with  the  predictions  of  the  theory  of  Southern 
and  Goodings,  and  they  permit  an  estimate  of  the  values  of  certain  of  the 
fundamental  single-ion  and  two-ion  magnetoelastic  coupling  constants  of  terbium. 

The  estimated  coupling  constants,  however,  exhibit  an  unexpected  temperature 
dependence  not  accounted  for  by  the  theory.  It  is  believed,  however,  that 
measurements  over  a  greater  temperature  interval,  using  considerably  larger 
applied  magnetic  field  values,  would  probably  lead  to  better  agreement  between 
the  experimental  results  and  the  theory  of  Southern  and  Goodings.  Such  an 

extension  of  the  work  reported  here  ,  ,  .  . , 

ported  nere  io  planned,  and  a  similar  experiment  has 

been  started  for  the  case  of  paramagnetic  dysprcsium. 

1.  Theory 


11 


In  this  section,  that  part  of  the  theory  of  Southern  and  Goodings1 - 
pertinent  to  the  present  discussion  is  outlined.  These  authors  utilised  both  single¬ 
ion  and  two-ion  magnetoelastic  interactions,  and  they  took  account  of  the 
quadratic  nature  of  the  conventionally  defined  finite  lattice-strain  tensor 
since  the  present  interest  lies  only  in  the  elastic  constants  c  and  c  ,  which 
are  measured  by  observing  the  velocities  of  longitudinally  polarised  elastic 
waves  propagated  along  the  crystallographic  a-,  b-,  or  c-axes,  the  magnetoelastic 
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interaction  Hamiltonian  used  by  Southern  and  Goodings  is  truncated  here  to  include 
the  effects  only  of  these  longitudinal  strains. 

The  single-ion  Hamiltonian  for  the  magnetoelastic  interaction,  subject  to 
the  truncation  mentioned  above,  is  given  by 


H 


me 


"  M20  l  E  'lj20(i)  "  M202  ?En'2j20(i)  “  M22  \  E}J22 (i) '  (HI-1) 


where  the  M-coeff icients  are  the  single-ion  magnetoelastic  coupling  constants 
defined  by  Southern  and  Goodings,  the  E-components  are  the  components  of  the  finite- 
strain  irreducible  tensor,  and  the  J-operators  are  the  second-rank  irreducible 
tensor  angular-momentum  operators  defined  in  a  previous  section  and  in  Appendix  I. 
The  index  i  is  a  label  for  each  ion  in  the  system.  The  two-ion  magnetoelastic 
Hamiltonian  is  given  by  the  following  expression: 


//TI  =  l  ll11  (i,j) 
me  . L  .  me  ' J 


(III-2a) 


V  +  Wft  - 


l 

T  J  .  *  J.] 
3  1  J 


(lll-2b) 


In  these  equations,  i  and  refer  to  a  particular  pair  of  ions,  and  the  D-  . 
coefficients  are  the  two-ion  coupling  constants.  The  M-  and  D-  coefficients  are 
determined  by  the  strain  dependence  of  the  magnotocrystalline  anisotropy  energy  and 
the  exchange  energy,  respectively. 

Southern  and  Goodings  recognized  that  in  the  paramagnetic  phase  all  the  second- 
degree  angular-momentum  operators  appearing  in  both  the  single-ion  and  two-ion 
magnetoelastic  Hamiltonians  given  above  would,  to  a  good  approximation,  contribute 
to  the  internal  energy  of  the  crystal  to  a  degree  proportional  to  the  square  of 
the  reduced  magnetization  and,  hence,  in  the  paramagnetic  phase,  to  the  square  of 

the  applied  magnetic  field.  Following  Southern  and  Goodings,  the  following  quantities 
can  be  defined: 


.  a  a .  l 

bl  =  N  M2J  JCJ-1/2) 


b2  =  N  M202j(J-1/2) 
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bY  -  N  ^22J(J-l/2) 


a“.  =  n  i  y  Da  . 

11  9  a  1 1  ■; 


2  j  llij 


dY  -  N  i  l  DY  . 

2.11 

3 

2 

,&  Jr  a 

d21  "  N  2  \  °21i j 


,01  J  r*  Ot 

d12  ~  N  2  |  °12ij 


,ot  J  r  a 

1  d22  ■  N  2  l  °22i j  " 

j 

where  N  is  the  number  of  ions  per  unit  volume.  The  following  expression  then 

corresponds  to  the  one  obtained  by  Southern  and  Goodings  for  the  magnetoelastic 
contributions  to  c, ,  and  c. 


"'ll 


'33  ' 


Ac33  '  +  2B)  P2  (cos  G)  +  C  +  D  ]  F  (T)  H'< 


( 1 1 1  —  3 ) 


AC11  “  1  (A  "  B)  p2(cos  9>  +  C  -  D  +  E  sin20  cos  2<f )  F  (T)  H' 


(III-4) 


'11 


in  the  above  equations,  F(T)  =  [  Xe(T)/MQ]2,  where  x0 (T)  is  the  anisotropic 
susceptibility  and  Mq  is  the  satuaration  magnetization.  The  angles,  0  and* , 
are  the  polar  angles  of  the  direction  of  the  magnetic  field  with  respect  to  the 
crystallographic  c-axis  and  a-axis.  In  eg.  (III-4)  the  plus  sign  is  used  when  c, 
is  to  be  determined  by  measuring  the  elast.ic-wave  velocity  along  the  a-axis. 

Thus,  according  to  the  calculations  of  Southern  and  Goodings,  the  changes  in 
the  two  elastic  constants  of  interest  in  an  applied  magnetic  field  should  be  propor¬ 
tional  to  the  square  of  the  value  of  the  applied  field  and  to  the  square  of  the 
magnetic  susceptibility.  The  functions  A  through  E,  defined  above,  all  functions 
of  the  fundamental  magnetoelastic  constants,  should  be  reasonably  constant  throughou* 
the  paramagnetic  temperature  range.  Examination  of  Eqs.  (III-3)  and  (III-4) 
reveals  that  if  the  magnetoelastic  changes  in  cn  and  c33  can  be  measured  accurately 
for  any  arbitrary  direction  of  the  applied  field,  then  the  constants  A  through 
E  may  all  be  determined  if  enough  independent  field  directions  are  chosen  for  the 
measurements.  The  actual  choice  of  directions  is  discussed  in  the  following  sec¬ 
tion  . 

2.  Experimental  Procedure 


The  single-crystal  specimen  of  terbium  used  in  this  work  was  prepared 
by  standard  spark-erosion  techniques,  as  described  elsewhere26.  For  the  measure- 
ment  of  c33,  20-MHz  longitudinally  polarized  elastic  waves  were  propagated  along  the 
crystallographic  c-axis,  and  for  the  measurement  of  c^  longitudinal  waves  of 
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the  same  frequency  were  propagated  along  the  b-axis.  The  elastic  constants  were 
measured  to  a  high  degree  of  accuracy  using  a  pulse-superposition  spectrometer 
described  elsewhere  to  measure  the  velocities  of  the  elastic  waves.  The 
temperature  of  the  specimen  was  maintained  at  any  desired  value  in  the  range  of 
interest,  230  K  to  300  K,  by  means  of  a  continuous-flow  cryostat  and  a  temperature 
control  system  of  commercial  design.  The  specimen  was  positioned  in  an  electro¬ 
magnet  such  that  the  field,  with  a  maximum  value  of  18  kOe,  could  be  directed 
along  any  desired  crystallographic  direction.  In  order  to  obtain  a  sufficient 
number  of  independent  relations,  employing  eqs.  (III-3)  and  (IH-4)  ,  to  permit 
a  complete  determination  of  all  the  quantities  A,  B,  C,  D,  and  E ,  three  different 
field  directions,  along  the  a-,  b-,  and  c-axes,  were  employed  in  the  determination 
of  the  field  dependences  of  the  two  elastic  constants  of  interest.  The  field 
orientation  used  in  the  measurement  of  each  elastic  constant  and  the  resulting 
relationaship  for  the  quantities  A  through  E  are  given  in  Table  I. 


TABLE  I 


Elastic  Constant 

Field  Direction 

Ac/F(T)H2 

r* 

33 

• 

c-axis 

A  +  2B  +  C  +  2D 

C33 

a-axis 

~j(A  +  2B)  +  C  +  2D 

C11 

c-axis 

A  -  B  +  C  -  D 

C11 

a-axis 

-  j(A  -B)  +  C-  D+  E 

C11 

b-axis 

-  j(A  -B)  +  C-  D-  E 

For  each  of  the  orientations  of  the  field,  the  field  dependence  of  the  elastic- 
wave  velocities  along  the  c-axis  (leading  to  the  determination  of  c33)  and  along 
the  b-axis  (leading  to  the  determination  of  cn)  were  measured  for  temperatures  from 
230  K  to  300  K.  Although  the  orientation  of  the  crystal  was  fairly  accurately 
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determined  through  the  use  of  X-ray  techniques,  the  measured  variation  of  the 
elastic-wave  velocity  with  respect  to  the  field  direction  could  be  used,  as  seen 

in  Fig.  17,  to  confirm  the  orientation  with  the  specimen  in  place  inside  the 
cryostat. 


Fig.  17  Angular  variation  of  b-axis  longitudinal-wave  velocity  as  field 
is  rotated  in  a-c  plane. 

3 .  Results 

The  results  of  this  work  are  depicted  in  Figs.  18  -  22,  where  the 
charges  in  the  longitudinal  elastic-wave  velocities  used  to  measure  the  elastic 
constants  of  interest  are  shown  as  functions  of  the  square  of  the  applied  magnetic 
field  for  various  temperatures.  The  straight  lines  shown  on  these  figures, 
tangent  to  the  actual  experimentally  determined  curves  at  low  fields,  were  used 
in  the  computations  discussed  below  to  represent  slope  (Av/v)/H^.  It  can  be 
seen  that  the  linear  dependence  of  the  relative  change  in  velocity  on  the 
square  of  the  applied  field,  as  predicted  by  the  theory  of  Southern  and  Goodings, 
is  in  reasonable  agreement  with  the  data  for  temperatures  well  above  the  Neel 
point.  At  the  lower  temperatures,  near  TN,  the  experimental  data  at  high  fields 
deviate  from  the  straight  lines,  undoubtedly  because  of  the  onset  of  saturation 
effects  in  the  magnetization. 
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Fig. 


18  Variation  of  c-axis 


longitudinal -wave  velocity  with  H  along  c-axis. 


Fig.  19  Variation  of  c-axis  longitudinal-wave  velocity  with  H  along  b-axis. 


The  results  shown  in  Figs.  18-22  were  analyzed  according  to  the  theory  of 
Southern  and  Goodings  in  the  nanner  described  above.  The  data  of  Hegland  et  ,1 

Worn  nt?A^  Pav  i-L _  _  .  •  <  .  .  .  2^  “  — 


32 


were  used  for  the  susceptibility,  and  the  data  of  Salama  et  al  26  for  the’ 
elastic  constants  in  zero  magnetic  field  were  employed  in  this  analysis,  s 
straightforward  numerical  computation  was  carried  out,  using  the  relations  of 
table  I  toobtain  values  for  the  quantities  »  through  E  at  each  temperature.  The  valu 
of  the  magnetoelastic  coupling  constants  were  then  computed  at  temperatures  from 
240  K  through  300  K.  and  the  results  are  presented  in  Table  II. 


DAAH01-72-C-0285  FINAL  REPORT 


42. 


Fig.  20  Variation  of  b-axis  longitudinal-wave  velocity  for  H  along  c-axis 

«*■  rtNMuM  Ti] 


Fig.  21  Variation  of  b-axis  longitudinal-wave  velocity  for  H  along  a-axis. 

It  can  be  seen  from  Table  II  that  there  is  a  substantial,  and  unexpected, 
variation  of  the  magnetoelastic  constants  with  temperature  over  the  range  240  K 
to  300  K.  This  variation  of  what  were  expected  to  be  nearly  temperature- 
independent  constants  over  this  limited  range  of  temperatures  is  illustrated  in 
Fig.  23,  where  the  quantity  3^/5  ♦  cT  is  omitted  for  the  sake  of  clarity. 

The  theory  of  both  the  single-ion  and  two-ion  magnetoelastic  interactions  in 
rare-earth  materials  is  based  on  the  assumption  that  the  crystal-field  and 
exchange  interactions  are  dependent  upon  the  lattice  strain7.  This  theory  has 
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TEMPERATURE  (K) 

Fig.  23  Temperature  dependence  of  the  magnetoelastic  coupling  constants  of 
single-crystal  terbium. 

been  quite  successful  in  explaining  many  of  the  observed  features  of  the  static 
magnetostriction  of  rare-earth  materials18.  One  result  of  the  comparison  of  the 
usual  theory  of  the  static  magnetoelastic  interactions  with  measurements  of 
static  magnetostriction  is  that  the  magnetoelastic  coupling  constants  are,  in  fact, 
essentially  constant  over  a  very  large  temperature  range,  extending  from  tempera¬ 
tures  in  the  liquid-helium  range  up  to  as  high  as  room  temperature.  It  is,  therefore, 
rather  surprising  that  such  a  large  variation  as  that  seen  in  Table  II  and  Fig.  23 
takes  place  over  the  limited  temperature  range  from  240  K  to  300  K.  One  would 
expect  only  a  very  small  variation  in  both  the  M-constants  and  the  D-constants  over 
this  range  and,  as  a  result,  only  a  very  small  variation  in  the  b-  and  d-quantities. 

The  theory  of  Southern  and  Goodings,  however,  upon  which  the  analysis  of  the 
data  obtained  in  the  present  work  was  based,  does  not  appear  to  have  neglected 
any  important  contributions  to  the  field  dependences  of  c^  and  c^  in  the  para¬ 
magnetic  state.  The  most  probable  reason  for  the  large  temperature  dependence  of 
the  results  given  in  Table  II  and  Fig.  23  is  the  fact  that  che  measurements 
reported  here  were  all  made  over  a  range  of  temperatures  fairly  close  to  the 
ordering  temperature,  Tn  =  229  K.  The  effects  of  short-range  order  near  a  phase 
transition  always  complicate  the  interpretation  of  the  observed  properties  of  any 
magnetic  material,  and  it  would  appear  that  such  is  the  case  here,  in  the  absence 
of  any  demonstrable  deficiency  in  the  theory  of  Southern  and  Goodings11. 
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It  would,  therefore,  be  desirable  to  attempt  to  make  similar  measurements  at 
somewhat  higher  temperatures,  where  the  effects  of  the  phase  transition  should  be 
less  important.  Such  measurements  at  higher  temperatures  would,  however,  be  diffi¬ 
cult  without  employing  much  larger  values  of  the  applied  field  than  were  available 
for  the  work  reported  here.  It  is  anticipated  that  this  work  will  be  extended 
to  higher  temperatures  in  fields  up  to  55  kOe,  and  through  such  additional  work 
it  may  be  possible  to  measure  the  magnetoelastic  coupling  constants  to  greater 
accuracy  than  possible  in  the  work  reported  here  and  to  verify  the  theoretical 
predictions  of  Southern  and  Goodings  in  a  more  satisfactory  manner. 

yit^sonic  Generation  in  Terbium-Iron  Thin  Films 


This  section  is  concerned  with  the  results  of  an  investigation  of  magneto- 
strictive  ultrasonic  generation  in  thin  films  of  terbium-iron  compounds,  particu- 

lHrly  CUbiC  Laves“Phase  compound  TbFe2,  which  was  recently  found,  by  Clark  and 
Belson  ,  to  exhibit  very  large  static  magnetostriction  at  room  temperature  in 
polycrystal line  specimens.  The  work  reported  here  was  motivated  by  earlier  studies 
on  magnetostrictive  ultrasonic  generation  in  rare-earth  materials24'25. 

In  the  cases  of  all  the  materials  studied  in  the  work  reported  here,  namely 
TbFe2,  TbFe^,  TbFe ,  and  Tb0 # 3°yo. Ve2'  filna  of  two-micrometer  approximate  thick¬ 
ness  exhibited  high  for  ultrasonic  generation  of  shear  elastic  waves  at  frequencies 
in  the  vicinity  of  680  MHz  at  room  temperature.  In  the  best  cases,  the  magneto¬ 
strictive  generation  was  about  100  times  more  efficient  than  that  for  piezoelectric 
longitudinal-wave  generation  in  single-crystal  X-cut  quartz.  A  very  interesting 
fe-ture  of  the  TbFe2  films  was  the  discovery  that  they  exhibited,  under  most  cases 
of  deposition  onto  quartz  substrates,  a  strong  uniaxial  anisotropy,  invariably 
with  an  easy  axis  of  magnetization  perpendicular  to  the  plane  of  the  film.  A 
large  remanent  magnetization,  amounting  to  some  80  per  cent  of  the  saturation 
value  of  the  magnetization,  was  often  observed,  and  coercive  forces  as  high  as 
5  kOe  were  measured.  As  a  result  cf  this  strong  anisotropy,  which  was  not  at  all 
expected,  in  view  of  the  rather  isotropic  nature  of  bulk  polycrystalline  TLFe  , 
strong  ultrasonic  generation  was  possible  even  in  the  absence  of  a  biassing  magnetic 
field.  An  interesting  aspect  of  the  existence  cf  the  anisotropy  is  the  result  that 
it  is  relatively  simple  to  prepare  a  thin-film  permanent  magnet  using  TbFe  ,  and  it 

is  possible  that  such  magnets  might  find  interesting  applications  in  areas  other  than 
those  considered  here. 
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1 •  Experimental  Technique 

The  specimen  films  were  rapidly  deposited  from  starting  bulk  material, 
prepared  by  ar-melting  in  an  argon  atmosphere,  onto  single-crystal  x-cut  quarts 
substrates  which  served  as  acoustic  delay  lines.  Deposition  was  carried  out  at 
substrate  temperatures  in  the  vicinity  of  300-  c  in  a  conventional  vacuum  system 
of  st-rting  pressure  approximately  lo’7  Torr.  The  evaporant  was  directly  heated  in  a 
tungsten  boat,  with  no  apparent  interaction  with  the  boat,  such  as  might  be  expected 
in  the  case  of  iron-group  materials.  The  film  thickness  for  which  most  measurements 
were  made  was  approximately  2,  pm,  as  determined  by  a  commercial  quartz-crystal 
thickness  monitor.  The  TbFe,  material  used  for  most  of  this  work  was  obtained  from 
a  commercial  source.  Some  TbFe2  and  all  the  remaining  materials  were  prepared  in 
this  laboratory,  with  the  exception  of  the  Tb-Dy-Fe  ternary  compound,  which  was 
kindly  provided  by  A.  E.  Clark  of  the  Naval  Ordnance  Laboratory. 

The  pulse  spectrometer  consisted  of  a  commercial  transmitter-receiver 
capable  of  operation  over  the  frequency  range  300-700  MHz,  and  it  was  equipped 
with  a  gated  "boxcar"  integrator  to  permit  the  detection  of  any  selected  ultrasonic 
echo  and  the  accurate  measurement  of  its  amplitude.  The  specimen  film  was  placed 
inside  a  stripline  cavity,  with  the  rf  magnetic  field  parallel  to  the  plane  of  the 
film.  A  dc  magnetic  field  as  large  as  18  kG  could  be  applied  using  a  conventional 
electromeagnet.  All  measurements  reported  here  were  made  at  room  temperature,  in 
the  vicinity  of  295  K.  The  magnetization  measurements,  described  below,  which  were 
carried  out  on  the  TbFe.,  films,  were  made  using  a  Foner  vibrating-sample  magneto¬ 
meter  of  high  sensitivity,  expecially  constructed  for  this  work. 

2.  Experimental  Result-.?; 


The  results  presented  here  on  magnetostrictive  ultrasonic  generation 
in  Tb-Fe  materials  were  obtained  at  an  ultrasonic  frequency  of  680-690  MHz  at 
room  temperature.  Ultrasonic  elastic  waves  were  generated  only  with  the  applied 
dc  magnetic  field  perpendicular  to  the  plane  of  the  film,  and,  in  all  cases,  the 
waves  were  transversely  polarized,  with  the  plane  of  polarization  parallel  to  the 
direction  of  the  rf  magnetic  field.  The  determination  of  tne  plane  of  polarization 
was  made  from  an  examination  of  the  shear-wave  modes  propagating  in  the  quartz 
delay-line  substrate.  The  substrate  was  normally  oriented  in  the  cavity  such  that 
the  faster  of  the  two  shear-wave  modes  was.  used  to  observe  the 


generation  and  detec- 
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of  the  elastic  waves.  In  the  figures  which  follow,  the  dc  magnetic  field  was  swept 

from  zero  up  to  a  maximum  value  of  18  kOe,  starting  in  the  direction  in  which  a  field 

was  last  previously  applied.  It  was  then  swept  back  through  zero  and  all  the  way  to 

16  koe  in  the  opposite  direction,  returning  finally  to  zero,  arrows  on  the  figures 
indicate  this  sequence. 


Fig. 24  Ultrasonic  echo  amplitude  as  a  function  of  applied  magnetic  field 
T^Fe2  on  quartz  substrate. 


The  field  dependence  of  magnetostrictive  ultrasonic  generation  in 
tbFe.,  films  is  shown  in  rig.  24.  This  material  was  selected  as  the  initial 
substance  to^be  studied  because  of  its  large  static  magnetostriction  at  room 
temperature  .  It  exhibits  ultrasonic  generation  efficiencies  much  higher  than 
that  of  any  other  magnetostrictive  materials  which  has  been  reported24'25.  The 
generation  efficiency  for  shear  elastic  waves  at  room  temperature  and  a  frequency 
of  680  MHz  was  typically  100  times  that  observed  for  longitudinal-wave  piezoelec¬ 
tric  generation  in  single-crystal  quartz  under  the  same  conditions.  The  large 
hysteresis  shown  in  Fig.  1  was,  however,  totally  unexpected,  since  it  implies 
a  coercive  force  of  the  order  of  4  kOe  for  these  thin  films.  Although  Clark 
ot  al  have  shown  that  TbFe2  does  possess  a  large  magnetic  anisotropy,  bulk 

specimens  of  this  material  do  not  exhibit  any  appreciable  hysteresis  nor  do  they 
possess  a  large  coercive  force. 

A  vibrating -sample  magnetometer  was  used  to  measure  the  magnetization  of  the 
2  films  used  in  this  work  and  to  compare  this  magnetization  with  that  of  the 
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bulk  material  utilized  for  the  preparation  of  the  films.  The  results  of  these  measure 
ments  are  shown  in  Figs.  25  and  26. 


Pig.  25  Magnetization  of  TbFe.,  Thin  films  vs  applied  magnetic  field. 


Pig.  26  Magnetization  of  bulk  polycrystalline  TbFe2  vs  applied  magnetic  field. 


The  magnetic  anisotropy  of  the  films  is  apparent  in  Pig.  25.  For  the  case  in 
which  the  field  is  applied  perpendicular  to  the  plane  of  the  film,  the  coercive 
force  is  approximately  5  kOe,  and  saturation  is  apparently  complete  at  a  field 
of  the  order  of  10-12  kOe.  For  the  case  in  which  the  applied  field  is  parallel  to 
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the  plane  of  the  film,  the  coercive  force  is  only  of  the  order  of  3  koe,  and 
saturation  is  not  reached  even  for  fields  of  the  order  of  20  koe.  By  contrast, 
the  specimen  of  bulk  TbFOj,  although  unsaturated  at  20  kOe,  exhibits  little  hys¬ 
teresis,  and  it  has  no  preferred  anisotropy  axis.  Clark34  has  reported  the 
development  of  high  coercive  forces  in  sputtered  amorphous  films  of  TbFe  after 
annealing  in  a  strong  magnetic  field  at  temperatures  above  300"  c.  The  source 
of  this  strong  anisotropy,  with  easy  axis  perpendicular  to  the  plane  of  the  film 
is  not  understood  at  this  time,  but  it  is  undoubtedly  due  to  at  least  some 
appreciable  degree  of  crystallographic  ordering.  The  anisotropy  dees,  however, 
yield  the  fortunate  result  that  a  TbFe.,  film  will  remain  strongly  magnetised  in 
the  absence  of  an  applied  magnetic  field,  after  first  being  magnetized  to 
saturation.  As  a  consequence,  it  is  capable  of  efficient  ultrasonic  generation  with 
no  bias  magnetic  field.  The  magnetization  data  of  Fig.  25  were  used  to  plot  the 
ultrasonic  generation  amplitude  of  Fig.  24  as  a  function  of  specimen  magnetization, 
rather  than  applied  field,  with  the  result  shown  in  Fig.  27.  Here,  it  can 
be  seen  that  the  ultiasonic  amplitude,  in  decibels,  is  nearly  linearly  dependent 
upon  the  specimen  magnetization,  a  result  predicted  to  some  extent  by  the  simple 
theoretical  model  of  ultrasonic  generation  presented  below. 


M  (#mu/gm) 


Fig.  27  Ultrasonic  echo  amplitude  as  a  function  of  film  magnetization;  TbFe., . 


0.3Dy0.7Fe2  was  also  investigated,  since  Clark  and 
material  also  exhibits  large  magnetostriction,  but 
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that  its  anisotropy  is  much  smaller  than  that  of  TbFe  • 

thin  fiw  «...  IDFe2*  u^trasonic  generation  in 

fU»s  this  material  is  shown  in  Flg.  28,  uhere  it  ca„  be  seen  ^ 

is  almost  no  hysteresis.  The  maxim™  efficiency  of  generation  in  this 

hown  xn  Fig.  .9.  Although  there  was  substantial  hysteresis  thpr* 

.  Hysteresis,  there  was  no  remanent 

ultrasonic  generation,  and  the  peak  eff f-  - 
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■Wave  Generation 


»  simplified  version  of  the  theoretical  model  of  Section  II  is  ,  , 

for  the  calculation  section  II  iS  employed  here 

-  h  materialsT  ^ 

implr  yr„g  pti°ns  discussed  ^  ^  ^  ^ 

sahe  avoiding  unnecessary  complexity,  that  the  rare-earth  ions  are  arranged 
simp  e  erromagnetic  array,  with  every  ion  equivalent  to  all  the  others 

t  present  purpose'  °niy  the  sin9ie-io-  ~ 
of  th  u  a“°“  °"ly  a  qUalitati''e  understanding  of  the  dependence 

1  “ltraSOniC  9e"eritl°"  »  —  magnetization  is  solght  at  this 

It  is  assumed  that  the  shear 

_  ‘  elastic  waves  propagate  along  the  z-axis  with 

polarization  in  the  x-v  n1a„e  in  .  y  2  axis  with 

y  P  e  m  an  isotropic  medium.  A  dc  maqnetir  f-ioia  • 

:  mr:;?::.::  i  rr- — - i! — - 

the  eW  •  he  Particle  displacement  in  the  medium  due  to 

the  elastic  wave  is  denoted  by  u(z,t)  for  th*  v 

y- component,  with  corresponding  shear  strain  e  ^  f°r  ^ 

singleton  magnetoelastic  Hamiltonian  per  ion  may  be  writ^5  =  9U/3Z'  ^  ^ 

Hme  ~  e4 ^JyJz  +  +  e  (J  j  +  j  j  )  ] 

y  2  y  5  x  z  zV  J  '  (III-5) 

where  B  is  a  magnetoelastic  coupling  constant,  and  the  J  are  the 
the  angular-momentum  operator  of  the  1  exponents  of 

tonian  neal  * ■  rare-earth  ion.  The  purely  magnetic  Han.il- 

an,  neglecting  anisotropy  energy  and  assuming  a  molecular- field  ■  . 

for  the  exchange  int  molecular- field  approximation 

the  exchange  interaction,  may  be  written  (per  ion) 


Hm  -  ^  VHrfJx  t  (H,  +  He)  ,y  , 


(III-6) 


^ere  n  is  the  time-dependent  appiied  rf  magnetic  field  along  the  x-axis  and 

°  “  1  6  aPPHed  dc  magnetic  field  along  the  z-axis.  The  field  II  is  the’ 

molecular  field,  which,  for  a  simrn  e  h  effective 

be  expressed  as  „  .  r  h  ,  "a5netized,  ferromagnet,  can 

molecular-field  constant  ^  'a  ^  Sat“rati°"  "“3"utization,  r  is  a 

constant,  and  c  is  the  reduced  magnetization,  M/M  .  To 

second  order  in  perturbation  theory  the  into  i  ° 

___  .  y'  th  lnternal  energy  associated  with  H 

can  be  expressed  in  the  following  way  me 
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U 

me 


B  e  <J  >H  ,/(H 
5  20  rf  o 


+  H  ) 
e 


(III-7) 


.  2 
In  this  expression,  the  operator  J2Q  is  defined  as  3J  -  J(J+1),  and  its 

thermal-equilibrium  expectation  value  is  given,  from  the  work  of  Callen  and 

Callen^,  as 


<J2Q>  =  J(2J-1)I5/2(W)/I1/2(W)  ,  (III-8) 

where  W  represents  the  inverse  Langevin  function  of  the  reduced  magnetization. 

Now,  if  U  is  regarded  as  a  strain-  and  field-dependent  contribution  to 
the  internal  energy,  to  be  added  to  the  normal  elastic  energy  density,  then 
the  following  wave  equation  for  shear  elastic  waves  propagating  along  the 
z_axis  and  polarized  along  the  x-axis  (parallel  to  the  rf  magnetic  field)  may 
be  derived  by  standard  techniques: 


p  3ju 
3t2 


=  c 


a2u 
44  3z2 


[NB<J  >/  (H  +  H  )  ]  -5—  H 

20  o  e  3z  rf 


(III-9) 


In  the  above  expression,  p  is  the  density  of  the  material,  and  c,.  is  a  normal 

elastic  constant.  The  second  term  on  the  right-hand  side  of  the  wave  equation 

represents  an  inhomogeneous  driving  term  in  the  wave  equation  which  is  responsi- 

bel  for  the  generation  of  the  elastic  waves.  The  z-dependence  of  arises,  of 

course,  from  the  finite  skin  depth  of  the  metallic  specimen.  Without  going  into 

the  details  of  the  solution  of  the  wave  equation,  it  can  be  seen  that  the 

amplitude  of  any  elastic  waves  excited  by  the  presence  of  the  driving  term  will 

be  proportional  to  the  bracketed  part  of  that  term.  Since  H  is  normally  much 

e 

greater  than  K^,  the  magnitude  of  the  driving  term  in  the  wave  equation  will  be 
proportional  to  <J2q>/  0  '  A  Plot  of  this  quantity  is  given  in  Fig.  30,  where 
it  can  be  seen  that  the  predicted  amplitude  of  the  magnetostrictively  generated 
elastic  waves  is  expected  to  vary  almost  linearly  with  the  specimen  magnetiza¬ 
tion.  Of  course,  the  fact  that  the  observed  generation  amplitude  shown  in 
Fig.  27  varies  linearly  as  the  specimen  magnetization  represents  a  somewhat  differ 
ent  effect,  being  due  to  the  rotation  of  domains  into  the  direction  of  the  applied 
field. 

It  should,  of  course,  be  recognized  that  the  present  treatment  is  simplified 
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to  the  point  at  which  much  detail  would  be  lost.  For  example,  any  possibility  of 
magnetic  resonance  has  been  excluded  from  the  formulation  given  here,  and  this 
would  seem  to  be  required  to  account  for  the  resonant-like  behavior  of  the  Tb-Dy-Fe 
compound  shown  in  Fig.  28  .  However,  the  experimental  result  that  the  polarization 
of  the  elastic  waves  generated  magnetostrictively  is  parallel  to  the  rf  field 
is  correctly  predicted  by  this  simplified  treatment,  and  it  might  be  expected 
that  a  fuller  treatment,  eliminating  s^me  of  the  large  number  of  simplifying 
assumptions  employed  here,  would  lead  to  a  reasonably  good  agreement  between 
experimental  results  and  theoretical  predictions. 

D.  Magnetoelastic  Effects  in  Erbium 

As  noted  on  numerous  occasions  in  the  preceding  sections  of  this 
report,  the  very  strong  magnetoelastic  interactions  found  in  the  heavy  rare- 
earth  elements  give  rise  to  a  number  of  interesting  phenomena.  The  subject  of  this 
section  is  a  description  of  two  dynamic  magnetoelastic  effects  observed  in  the  course 
of  studies  of  the  general  characteristics  of  elastic-wave  propagation  in 
single-crystal  erbium.  Erbium  is  the  first  of  the  rare-earth  series  of  elements 
all  of  whose  ordered  magnetic  phases  include  ordered  c-axis  components  of  the 
magnetization.  For  this  reason,  it  is  possible  to  observe  the  phenomenon  of  cir¬ 
cular  magnetoacoustic  birefringence  in  erbium,  and  this  section  reports  on  the 
observation  of  this  effect.  It  also  occurs  that  the  dependence  of  longitudinal 
dlastic-wave  velocities  on  the  value  of  an  applied  magnetic  field  is  quite  large, 
and  this  section  reports  on  this  effect,  manifested  as  a  very  large  dependence  of 
the  elastic  constant  c^  on  the  strenght  of  the  applied  field. 

Specifically,  the  work  described  in  this  section  concerns  the  following  two 

magnetoelastic  effects.  The  first  is  the  anomalous  temperature  dependence  of  rhe 

elastic  constant  c.^,  which  undergoes  a  step-like  decrease  at  the  magnetic  phase 

26—30 

transitions  rather  than  the  cusp-like  behavior  exhibited  by  other  elements 
Furthermore,  the  application  of  a  magnetic  field  of  30-50  kG  produces  a  very  large 
decrease  ii  the  value  of  c33  in  the  vicinity  of  T^  (86  K) .  This  large  decrease, 
amounting  to  approximately  9  per  cent,  is  much  larger  than  any  of  the  zero-field 
elastic-constant  anomalies  which  have  been  reported.  The  second  effect  reported 
here  i.s  the  observation  of  a  large  circular  birefringence  for  elastic  shear  waves 
propagating  along  the  c-axis  of  Er.  At  a  temperature  of  40  K  in  an  applied  c-axis 
field  of  15  kG,  this  effect  leads  to  a  rotation  of  the  plane  of  polarization  of 
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approximately  0.65  rad/ cm  at  a  frequency  of  20  MHz.  An  attempt  has  been  made  to 
explain  the  results  given  here  in  terms  of  the  theoretical  model  of  Section  II 
of  this  report.  The  complicated  magnetic  structures  of  erbium  make  detailed  cal¬ 
culation  difficult,  but  some  degree  of  qualitative  agreement  between  theory  and 
experiment  has  been  achieved. 

1.  Experimental  Technique 

The  single-crystal  erbium  specimens  used  in  this  work  were  prepa  ‘ed  by 
standard  spark-erosion  techniques,  with  flat  and  parallel  end  surfaces  perpen¬ 
dicular  to  the  c-axis.  The  length  along  the  c-axis  was  approximately  '  m.  For  the 
measurement  of  c^»  the  pulse-superposition  method  was  used  to  deter.  <.  the 
elastic-wave  velocity  to  high  precision.  For  the  birefringence  measurements,  a 
Matec,  Inc.  pulsed  transmitter-receiver  was  employed  to  permit  the  measurement  of 
the  amplitude  of  a  single  ultrasonic  echo.  The  ultrasonic  frequency  was  20  MHz 
for  both  types  of  measurement.  For  the  measurement  of  c33  longitudinally  polarized 
waves  were  employed,  whereas  for  the  observation  of  the  birefringence  transversely 
polarized  waves  were  used.  In  the  case  of  circular  birefringence,  which  leads  to 
a  rotation  of  the  plane  of  polarization  of  a  linearly  transversely  polarized  wave 
as  it  propagates,  this  rotation  led  to  an  oscillatory  dependence  of  the  detected 
echo  amplitude  as  the  applied  magnetic  field  was  varied,  causing  the  rotation 
per  unit  length  to  change.  An  example  of  this  oscillatory  behavior  is  shown  in 
fig.  31.  For  the  measurement  of  c33  and  its  dependence  on  applied  fields,  a 
superconducting  solenoid  capable  of  fields  up  to  60  kG  was  employed.  (For  the 
birefringence  measurements,  a  conventional  iron  electromagnet  capable  of  fields  up 
to  18  kG  was  used.)  For  all  measurements  the  specimen  temperature  was  controlled 
within  .05  K  of  the  desired  temperature  through  the  use  of  a  continuous— flow 
cryostat  and  conventional  sensor-heater-feedback  control  system. 

2.  Circular  Magnetoacoustic  Birefringence 

For  the  observation  of  circular  birefringence,  it  is  necessary  that  the 
elastic  waves  propagate  along  a  crystallographic  direction  for  which  the  two 
shear  modes  are  normally  degenerate  in  the  absence  of  the  magnetoelastic  interaction 
It  is  also  necessary  that  there  be  a  component  of  the  magnetization,  either 
spontaneous  or  induced,  along  the  axis  of  propagation.  Erbium  is  the  only 
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readily  available  rare-earth  element  for  which  these  requirements  can  be 
conveniently  met,  in  that  all  of  its  magnetic  ordering  phases  (discussed  in 
Section  II,  Part  A)  involve  ordered  magnetization  components  along  the  c-axis. 
The  c-axis  is,  of  course,  the  only  axis  in  a  hexagonal  crystal  for  which  the  two 
shear  modes  are  degenerate.  Figures  31  and  32  show  some  of  the  results  obtained 
in  the  attempt  to  measure  circular  birefringence  in  single-crystal  erbium. 


Fig.  31  Oscillations  in  the  amplitude  of  the  n—  shear-wave  echo  propagating 
along  the  c-axis  of  single-crystal  Er. 


Fig.  32  Rotation  of  the  plane  of  polarization  of  shear  waves  propagating  along 
the  c-axis  in  single-crystal  erbium. 
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Fig.  31  shows  the  oscillatory  behavior  of  the  echo  amplitude  as  a  function  of 
the  applied  field  for  several  temperatures.  The  maxima  of  the  oscillations  occur 
when  the  rotation  of  the  plane  of  polarization  is  an  integral  multiple  of  it, 
and  the  minima  occur  for  odd  half-integral  multiples  of  it.  From  the  data  of 
Fig.  31,  the  curves  of  Fig.  32  were  derived,  giving  the  rotation  per  unit  length 
of  travel  along  the  c-axis  as  a  function  of  applied  field  for  several  temperatu-es 
in  the  so-called  quasi  antiphase  domain  region  of  erbium. 

The  observed  rotation  of  the  plane  of  polarization  for  shear  waves  in  erbium 
could  only  be  measured  for  temperatures  in  this  ordering  range  because  fields 
only  as  high  as  18  kG  were  available  for  these  measurements.  Higher  fields  are 
necessary  to  produce  a  substantial  c-axis  net  magnetization  component  at  higher 
temperatures.  As  a  result,  the  analysis  of  the  results  is  somewhat  more  complicated 
than  would  be  the  case  of  measurements  in,  for  example,  the  paramagnetic 
temperature  range.  Nevertheless,  the  theoretical  approach  described  in  Section  II 
has  been  employed,  again  in  a  very  simplified  form,  for  the  interpretation  of  the 
birefringence  results  in  terms  of  the  fundamental  magnetoelastic  interactions. 

An  estimate  for  the  shear  single-ion  magnetoelastic  coupling  constant  bZ'2 
of  100, K  per  ion  has  been  made,  and  this  value  is  somewhat  larger  than  the  value 
for  B.,'  of  20  K  per  ion  which  can  be  estimated  from  the  static  magnetostriction 
results  of  Rhyne  and  Legvold  5.  An  improved  agreement  between  theory  and  experi¬ 
ment  can  probably  be  obtained  when  measurements  at  higher  temperatures  are  carried 
out. 

t 

3*  Temperature  Dependence  of  c  for  Erbium 

The  temperature  dependence  of  c^  both  in  the  absence  of  an  appl^d 

magnetic  field  and  in  a  field  of  30  kG  along  the  c-axis  is  shown  in  Fig.  33. 

The  principal  behavior  in  zero  field  is  the  step-like  behavior  at  the  Neel 

point,  86  K  and  at  the  lower  phase  transition  at  52  K.  The  magnitude  of  the  step 

at  Tn  is  well  understood,  as  discussed  in  Section  III,  Part  A,  where  good  agreement 

between  theory  and  experiment  was  found  when  the  magnitude  of  this 

step  was  calculated.  The  application  of  a  field  less  than  that  required  to  induce 

a-  transition  from  the  sinusoidally  modulated  c-axis  oraeied  magnetization  to 

ferromagnetic  alignment  has  little  effect  on  c^,  but  above  the  threshold  for 

this  field-induced  AFM-FM  transition,  the  field  gives  rise  to  a  very  large  decrease 

in  the  value  of  c^.  Because  the  state  of  the  magnetization  below  T  is  not 

N 
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very  well  understood,  it  is  difficult  in  this  case  to  calculate  the  magnitude  of 
the  change  in  c33  due  to  the  presence  of  the  field.  If,  however,  the  ferromagnetic 
alignment  is  assumed  to  be  complete  just  below  in  the  presence  of  the  large 
applied  field,  the  theory  discussed  above  leads  to  a  predicted  decrease  in  c 
of  12  per  cent,  only  slightly  larger  than  the  observed  decrease,  if  the  magneto- 
elastic^  single-ion  coupling  constants  deduced  from  the  results  of  Rhyne  and 
Legvold  5  are  used  in  the  theory,  The  other  elastic  constants  of  erbium,  which 
also  exhibit  significant  zero-field  anomalies,  are  only  slightly  affected  by  the 

application  of  magnetic  fields  up  to  55  kG  either  parallel  to  or  perpendicular  to 
the  c-axis. 

E‘  Polo9raphic  Interferometric  Measurf-j-ent  of  Magnetostriction  in 
Rare-Earth-Iron  Compounds 

This  section  reports  the  results  of  the  application  of  the  relatively 
new  technique  of  holographic  interferometry  to  the  measurement  of  the  static 
magnetostriction  of  small  single-crystal  specimens  of  the  very  interesting 
cubic  Laves-phase  rare-earth-iron  compounds  such  as  TbFe2 .  As  previously  mentioned, 
TbFe2  exhibits  very  large  magnetostriction  in  polycrystalline  form  at  room 
temperature  ,  and  it  also  possesses  a  very  large  magnetic  anisotropy36,  although 
no  single-crystal  anisotropy  measurements  have  been  reported.  In  the  research 
program  described  here,  we  have  been  successful  in  growing  small  single-crystal 
specimens  of  two  RFe.,  materials,  namely  ErFe2  and  a  ternary  compound, 

Er0.314Tb0.686FV  Althou9h  these  crystals  were  too  small  to  permit  the  usual 
strain-gauge  type  of  measurement  of  the  magnetostriction,  it  was  possible, 
using  the  technique  of  holographic  interferometry,  to  measure  the  magnetostriction 
along  various  crystallographic  directions  for  these  materials,  and,  in  conjunction, 
measurements  of  the  magnetization  were  also  made. 

1.  Crystal  Growth 

The  crystals  employed  in  this  work  were  grown  by  electron-beam 
zone  refining  of  polycrystalline  material  which  had  been  previously  arc  melted 
from  the  appropriate  proportions  of  high-purity  specimens  of  the  constituent 
components.  In  general,  it  was  not  possible  to  obtain  single-crystal  specimens 
with  any  degree  of  repeatability,  and  the  results  presented  here  were  obtained 
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with  what  were  essentially  randomly  and  accidentally  grown  small  single-crystal 
grains  extracted  from  a  larger  polycrystalline  mass.  Because  of  the  peritectic 
behavior  of  most  rare— earth— iron  compounds  ,  it  will  undoubtedly  be  difficult  to 
obtain  large  single -crystal  specimens  grown  under  tightly  controlled  conditions, 
and  so  it  was  felt  that  it  would  be  worthwhile  to  carry  out  several  types  of 
experimental  mecisurement  before  attempting  to  grow  better  crystals.  The  details 
of  the  crystal-growing  procedure  are  given  elsewhere37. 

Because  the  compound  ErFe2  appeared  to  be  most  appropriate  for  the  growth 
of  single  crystals,  this  material  was  chosen  initially  as  the  subject  of  this 
investigation.  However,  the  first  lot  of  supposedly  pure  polycrystalline 
erbium  which  was  used  in  the  preparation  of  the  polycrystalline  ErFe2  actually 
contained  a  large  percentage  of  terbium.  As  a  result,  the  small  single-crystal 
specimen  which  was  obtained  from  this  material  was,  in  fact,  a  Tb-Er-Fe 
compound  having,  apparently,  the  same  cubic  Laves-phase  structure  as  TbFe2  or 
ErFe2.  The  composition  of  this  specimen  was  determined  by  the  technique  of 
electron-stimulated  X-ray  fluorescence,  using  an  electron  microprobe  apparatus. 

The  composition  of  the  crystal  actually  used  for  magnetostriction  and  magnetization 

measurements  was  Er0.3i4Tbo.686Fe2 ’  The  startin9  Er-Tb  mixture  did  not  contain 

nearly  as  high  a  proportion  of  Tb  as  did  the  single-crystal  specimen ,  but  it 

may  be  that  the  observed  Tb-Er  ratio  in  the  final  crystal  is  the  result  of 

some  degree  of  selectivity  in  the  actual  crystal-growth  process.  A  second 

lot  of  pure  erbium  was  used  in  the  preparation  of  a  second  single-crystal  specimen, 

and  its  X-ray  fluorescence  analysis  indicated  that  it  was  essentially  pure 

ErFe2 ■  crystals  were  irregularly  shaped,  with  linear  dimensions  ranging  from 

one  to  four  millimeters. 

2 .  Experimental  Technique  for  Holographic  Measurement  of  Magnetostriction 

Because  it  had  been  anticipated  that  only  very  small  crystals  of 
RF®2  comI?ounds  would  be  available,  making  the  use  of  conventional  methods  for  the 
measurement  of  magnetostriction  difficult,  if  not  impossible,  a  system  was  developed 
to  permit  the  measurement  of  small  deformations  of  small,  irregularly  shaped, 
bodies  using  the  method  of  real-time  holographic  interferometry.  Briefly,  the 

O  *7 

method,  which  is  described  in  considerable  detail  elsewhere  ,  may  be  described 
as  follows.  A  hologram  is  ofrmed  of  the  object  whose  deformation  is  to  be 
investigated,  using  coherent  light  from  and  argon- ion  laser.  The  hologram  is  then 
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subjected  to  the  same  light  source  as  that  used  in  its  original  formation  in 
such  a  way  that  the  holographically  reconstructed  three-dimensional  virtual  image 
of  the  object  coincides  with  the  object  itself.  Now,  when  the  object  is  displaced 
or  deformed  in  any  way,  light  reflected  from  its  surface  witll  set  up  an  inter¬ 
ference  pattern  with  light  from  the  holographically  reconstructed  image  of  the 
original  undeformed  object.  By  careful  measurements  on  the  positions  of  the  fringes 
which  can  be  seen  as  a  result  of  the  interference  pattern,  it  is  possible  to 
analyze  all  components  of  both  the  displacement  and  the  strain  of  the  object. 

In  the  present  application,  the  object  to  be  studied  with  holographic 
interferometry  is  the  magnetostrictive  RFe^  sin9le-crystal  specimen,  which  is 
located  between  the  poles  of  an  electromagnet.  A  hologram  is  first  formed  in  the 
absence  of  an  applied  magnetic  field,  and  its  reconstructed  image  is  superimposed 
upon  the  object,  resulting  in  a  fringeless  interference  pattern.  The  magnetic 
field  is  then  increased,  giving  rise  to  magnetostrictive  deformation  of  the 
RFe2  crystal  and  producing  a  pattern  of  light  and  dark  interference  fringes. 

Analysis  of  the  positions  of  these  fringes  then  permits  a  complete  determination 
of  all  the  modes  of  magnetostrictive  strain  in  the  specimen.  The  details  of 

37 

the  analysis  are  quite  complicated,  and  they  are  given  in  great  detail  elsewhere 

Measurements  were  also  made  on  the  specimen  magnetization  using  a  conventional 
For.er  '’ibrating-sample  magnetometer  constructed  in  this  laboratory 

3.  Experimental  Results 

The  static  room-temperature  magnetostriction  of  the  two  single-crystal 

specimens  along  various  crystallographic  axes  are  depicted  on  Figs.  34  and  35.  The 

magnetization  data  are  shown  in  Figs.  36  and  37.  Although  these  data  are  only 

preliminary,  they  reveal  quite  interesting  features.  Consider  first  the  case  of 

ErFe_.  The  magnetization  curves  are  quite  similar  to  those  given  by  Clark  and 
^  36 

Belson  ,  as  expected,  showing  a  high  degree  of  saturation  for  magnetization  along 
the  [111]  direction  both  in  the  case  of  the  magnetization  and  the  magnetostriction. 
The  magnetostriction  is  not  very  large,  amounting  to  at  most  a  few  hundred  parts 
per  million,  although  even  this  valur  is  considerably  larger  than  that  of  most 
iron-group  materials.  For  the  Er-Tb-Fe  ternary  crystal,  again  the  [111]  direction  is 
the  easy  axis,  but  the  anisotropy  is  apparently  much  stronger.  What  is  really 
interesting  is  the  magnitude  of  the  magnetostriction,  which,  for  the  [111]  direction 
has  a  saturation  value  of  X  =  5000  x  10  ^  as  large  as  that  of  many  pure  rare 
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earths  at  low  temperature.  This  value  is  also  much  greater  than  that  of  polycrysta- 

21  6 

lline  TbFe2  as  reported  by  Clark  and  Belson  (approximately  2000  x  10  ) .  Since 

terbium  and  erbium  are  expected  to  contribute  to  the  total  magnetostriction  with 
opposite  signs,  this  large  value  for  the  ternary  crystal  is  surprising. 

These  results  are  quite  preliminary,  and  they  have  not  been  fully  analyzed. 
They  are  presented  here  as  an  example  of  interesting  and  new  results  which 
will  probably  not  be  extended  or  further  analyzed  because  of  total  lack  of  funding 
for  this  entire  line  of  research. 


Fig.  33  Temperature  dependence  of  the  elastic  constant  c33  of  -single-crystal 
erbium  in  zero  applied  field  and  in  a  30-kG  c-axis  field. 


Fig.  34  Magnetostriction  of  Single-Crystal  ErFe2  at  room  temperature 
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Fig.  35  Magnetostriction  of  single-crystal  Er  Tb  Fe 

0.314  0.686  2 

room  temperature. 


Fig.  36  Magnetization  of  ErFe^  at  room  temperature. 


at 
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Fig.  37  Room  temperature  magnetization  of  Er  Tb  Fe 

0.314  0.686  2’ 


IV.  CONCLUSION 

This  report  has  presented  a  number  of  aspects  of  the  dynamical  magnetoelastic 
properties  of  rare-earth  materials.  It  will  be  noted  by  the  careful  reader  that 
in  most  cases  the  work  is  presented  in  a  preliminary  or  incomplete  state.  This 
situation  is  somewhat  unfortunate,  but  it  represents  the  state  of  affairs  a.t  the 
time  support  for  this  program  was  terminated.  In  this  program  it  was  necessary 
to  carry  out  a  great  amount  of  preliminary  work  before  results  began  to  come  in. 
During  the  last  year  of  the  project,  as  can  be  seen,  there  were  numerous  projects 

in  operation,  all  of  which  were  left  with  much  work  remaining  to  be  finished  when 

funds  for  the  support  of  this  work  were  terminated.  It  will,  of  course,  be  possi¬ 
ble  to  carry  on  with  the  analysis  of  many  of  these  results  without  financial 
support,  but  it  will  be  quite  some  time  before  all  these  results  are  completely 

analyzed  and  a  final  attempt  is  made  to  arrive  at  a  complete  and  coherent  picture 

of  dynamic  magnetoelastic  interactions  in  rare-earth  materials. 
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APPENDIX  I  DEFINITION  OF  IRREDUCIBLE-TENSOR  ANGULAR-MOMENTUM  OPERATORS 

Since  the  original  introduction  by  Stevens  of  the  so-called  operator 
equivalents  for  dealing  with  problems  of  crystal  fields  interacting  with  magnetic 
ions  in  crystal  lattices,  there  have  been  several  conventions  for  the  definition 
and  normalization  of  these  operators,  with  a  certain  degree  of  resulting  con¬ 
fusion  in  the  literature.  What  is  proposed  here  is  a  simple  definition  of  these 
operators  and  a  definition  which  is  not  only  easy  to  remember  but  which  makes  the 

computation  of  matrix  elements  of  these  operators  quite  simple.  The  definitions  which 

38 

follow  were  inspired  by  a  book  of  Watanabe 

We  shall  define  a  set  of  irreducible-tensor  angular-momentum  operators 

of  rank  i  as  J™  ,  where  m  takes  on  integral  values  from  -£  to  +Jt.  In  order  to 

arrive  at  a  consistent  set  of  operators  which  can  be  easily  derived  at  any  time, 

l 

we  define  the  operator  J  in  the  following  way: 


=  (-1) V  ;  J+  =  Jx  +  ijy  (AI-1) 

We  then  find  the  remaining  operators  of  this  rank  by  using  the  basic  commutation 
relations  obeyed  by  any  set  of  irreducible-tensor  operators: 

[  J.  ,  1  =  /  £U+1)  -  m(m-l)  (AI-2) 

With  this  definition  convention,  all  angular-momentum  operators  can  readily 
be  derived  and  their  matrix  elements  can  easily  be  evaluated.  As  an  example  of 
the  application  of  the  convention,  consider: 

J1  =  "J+  ;  J1  =  (2)1/2jz  ;  J11  =  J- 

J2  =  J+  ;  J2  =  "V2Jz  *  15  ;  J2  =  <2/3>1/2<3J2  "  J<J+U> 

j"2  *=  J2  ;  J"1-  J  (2J  -  1) 

z  —  z  —  z 

The  real  beauty  of  this  definition  becomes  apparent  when  numerical  computations 
must  be  made,  since,  if  a  matrix-oriented  computer  language  such  as  APL  is  available, 
the  matrices  for  all  desired  irreducible-tensor  operators  can  be  obtained  with 


l 
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great  ease.  The  principal  advantage  in  using  the  definition  given  here  is  that 
it  can  be  remembered  easily,  and  with  its  use  there  should  never  be  any  uncertainty 
about  the  definition  or  normalization  of  these  operators. 

In  addition  to  the  basic  irreducible  tensor  operators,  which  are,  of  course, 
given  for  spherical  irreducible  representations,  it  is  often  necessary  to  use 
operators  which  belong  to  irreducible  representations  of  lower  symmetry.  In  most 
cases  the  following  set  of  Hermitian  operators  will  be  useful: 


jL  -  ,  j-m  =  n-i>v:  -  ,  j,„  =  j? 

Examples  are: 


iO 


Jn  ■  J*  1  Ju  -  Jy  '  J10  =  JZ  12)1/2 


J22  =  JX  “  Jy  ;  J21  =  IJXJZ  +  Vx1  ;  etC* 


APPENDIX  II  IRREDUCIBLE-TENSOR  STRAIN  COMPONENTS  FOR  HEXAGONAL  SYMMETRY 


The  simplest  choice  for  the  second-rank  irreducible-tensor  strain 
components  for  hexagonal  symmetry  would  appear  to  be  the  following: 


e  + 


XX 


e  + 

yy 


e 

zz 


e 


a, 2 


3e 


zz 


-  e 


a,0 


e  -  e 
xx  yy 


eY2  =  e 


+  e 
xy  yx 


e  + 
xz 


e 


zx 


e  +  e 
yz  zy 


Using  other  definitions  with  various  multiplicative  factors, both  rational  and 
irrational,  which  stem  to  vary  from  one  author  to  another,  seems  to  cause 
confusion  and  makes  it  much  more  difficult  to  follow  through  calculations 
made  by  different  authors.  It  would  appear  to  be  to  the  advantage  of  all  new 
workers  entering  fields  where  these  irreducible-tensor  components  must  be  used 
frequently  to  standardize  the  definition,  and,  in  the  interests  of  simplicity, 
to  make  such  a  standard  definition  as  free  from  confusing  factors  as  possible. 
Some  authors  may  argus  that  the  multiplicative  factors  which  they  happened  to  use 
in  their  definitions  of  these  tensor  quantities  arise  naturally,  but  there  seems 
to  be  no  evidence  to  support  such  a  contention. 
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